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Inertial Navigation System Using Three 
TDF Gyroscopic Sensors Not Jointly Mounted 
On a Stable Platform 

S ummary 

An inertial navigation system is descirbed and analyzed based on two two- 
degree-of-fr eedom Schuler-gyropendulums and one two-degree-of-freedom 
azimuth gyro. The three sensors, each base motion isolated about its two 
input axes, are mounted on a common base, strapped down to the vehicle. 

The up and down pointing spin vectors of the two properly tuned gyro- 
pendulums track the vertical and indicate physically their velocity with 
respect to intertial space. The spin vector of the azimuth gyro is pointing 
northerly parallel to the earth axis. 

The system can be made self-aligning on a stationary base. If external 
measurements for the north direction and the vertical are available, 
initial disturbance torques can be measured and easily biased out. 

The error analysis shows that the system is practicable with today's 
technology, 

Tragheitsnavigationssystem unter Verwendung von drei, 
nicht auf einer gemeinsamen stabilen Plattform 
montierten Kreiselsensoren 

tibersicht 

Ein Tragheitsnavigationssystem, aufbauend auf zwei Schuler-Kreiselpendeln 
und einem Azimutkreisel mit je zwei Mefifreiheitsgraden, wird beschrieben 
und analysiert. Die drei Sensoren sind auf einer gemeinsamen Grundplatte 
montiert, die fest mit dem Fahrzeug verbunden ist, Dabei ist jeder Sensor 
fur sich urn seine beiden Eingangsachsen gegeniiber der Fahrzeugbewegung 
isoliert. 

Die auf- und abweisenden Drallvektoren der beiden speziell abgestimmten 
Kreiselpendel zeigen die Fahrzeuggeschwindigkeit gegeniiber dem Inertial- 
raum an und ermbglichen es, die Lotrichtung zu ermitteln, Der Drallvektor 
des Azimutkreisels weist nach Norden, parallel zur Erdachse, 

Eine vereinfachte Fehleranalyse ermbglicht, die Ausfuhrbarkeit des Systems 
zu beurteilen. 

Das System ist fahig zur Selbstausrichtung, Wenn externe Referenzen fur 
die Nordrichtung und die Vertikale vorhanden sind, kbnnen Anfangsstbr- 
momente gemessen und recht einfach kompensiert werden. 




4 




Acknowledgements 

This study was accomplished while the author was supported as a Postdoctoral 
Resident Research Associate by the National Academy of Sciences. Diiring his 
stay at NASA/Electronics Research Center the author was on a leave of 
absence from Deutsche Forschungs- und Versuchsanstalt fur Luft- und Raum- 
fahrt, Braunschweig, Germany. 

Dr. Arthur H. Lipton, the author's Scientific Advisor at NASA/ERC, gave 
his encouragement for this study, provided the atmosphere in which it could 
be carried out and gave valuable cirticism. Dr. Shigeo Okubo showed great 
interest in this study. The author likes to remember the many discussions 
he had with him. Professor Walter Wrigley, the author's honorable teacher 
at Massachusetts Institute of Technology, was open to consultation and 
provided valuable hints. The final analysis was reviewed by Dr. John E. Bortz, 
(NASA/ERC), Dipl.-Ing. Dietrich Rahlfs and Ing. Egmar Liibeck (both DFVLR) , 
from whom the author received estimable contributions. 

The author wishes to express his appreciation to all of them.; 

The publication of this report does not cinstitute approval of National 
Academy of Sciences, National Aeronautics and Space Administration or 
Deutsche Forschungs- und Versuchsanstalt fur Luft- und Raumfahrt of the 
findings or conclusions contained therein. It is published only for the 
exchange and stimulation of ideas. 


5 



TABLE OF CONTENTS . Page 

No. 

Abstract . ^ 

Acknowledgements 5 

List of Figures 8 

List of Tables 10 

1 . Notation .11 

2. Introduction 17 

3. Basic Performance and Implementation of the Navigation 20 

System 

4. Performance of the Gyropendulum in a Dynamic Environ- 
ment under Ideal Tuning and Compensation Conditions 31 

5. Performance of the Azimuth Gyro 44 

6. Design Considerations of the Control Loops for the 

Gyropendulum and the Azimuth Gyro 46 

6.1 Tuning-Control of the Gyropendulum 46 

6.2 Gimbal Control and Command Torques 51 

7. Simplified Error Analysis 54 

7.1 Computational Velocity Errors of a Gyropendulum 

with Signal Generators Aligned with the Naviga- 
tional Frame 55 

7.2 Attitude Errors and Reqiiirements for the Azimuth 

Gyro 60 

7.3 Frequencies of the Error Oscillations 63 

8. Initial Alignment Considerations 64 

8.1 Initial Alignment Based on External Measurements 65 

8.2 Self Alignment 68 

9 . Summary ^ 72 


6 


Appendix A Page 

No. 

Derivation of the Performance Equations for a Gyropendulum with 
Vertical Spin Axis in a Dynamic Environment 75 

A1 Readout Angles in the Navigational Frame 75 

A2 Readout /ingles in the Body Frame 83 

A3 Application of the Command Torques in the Free Azimuth Frame 84 

Appendix B 

Derivation of the Performance Equations of a Two-Degree-of- 
Freedom Azimuth Gyro with its Spin Vector Pointing Northerly, 

Parallel to the Earth Axis 86 

B1 Readout Angles in the Equatorial Tangent Frame 86 

B2 Readout Angles in the Body Frame for a Gimballed Gyro 87 

I 

B3 Readout Angles in the Body Frame for an Electrostatic 

Gyro (ESG) with Direction Cosine Pattern Readout 89 

Appendix C 

Derivation of the Error Equation of the Gyropendulum within the 
Navigational Frame 91 

References 99 


7 



f 


LIST OF FIGURES 

Figure Page 

No. No. 

3.1 A Gas Film Supported Free Rotor Gyro 22 

3.2 Gimballed Gyropendulum 23 

3.3 Gimballed Azimuth Gyro 26 

3.^ Signal Flow Diagram for the Computation of Azimuth I 4) ) 
and Latitude (L), Based on the Readout of the Azimuth 
Gyro and Attitude, Derived from the Gyropendulums 2 28 

3.5 Signal Flow Diagram for the Simplified Computation of 

Attitude (tp , ■&), Ground Speed and Position (L,l), 

Based on the ^uasi-Static Readout of 2 Schuler-Tuned 
Gyropendulums 29 

4.1 Signal Flow Diagram for a Two-Degree-of-Freedom 

Gyropendulum with Vertical Spin Axis 32 , 

4.2 Coning Motion of a Gyropendulum with Spin Vector Pointing 

Up or Down, Respectively 37 

4.3 Signal Flow Diagram for the Exact Computation of Attitude 

» Ground Speed (Vj^, V^) and Position (L,l), Based 
on the Readout of 2 Properly Tuned Gyropendulums 43 

6.1 Functional Diagram of a Navigational System Using Two 
Two-Degree-of-Freedom Gyropendulums and On* TDF 

Azimuth Gyro 49 

6.2 Signal Flow Diagram for the Gimballed Gyropendulum 50 

6.3 Angles between Coordinate Frames 51 

8.1 Signal Flow Diagram for Measuring 0^^^ and 0^,^ of a 

Gyropendulum During Initial Alignment Based on External 
Measurements 66 


8 


LIST OF FIGURES (CONTINUED) 

Figure Page 

No» No. 

A1 Model of the Pendulous Two-Degree of Freedom Gyroscope . 

with Vertical Spin Axis 75 

A2 The Navigational Frame and the Gyro Frame 

B1 The Equatorial Tangent Coordinate Frame S6 

B2 Readout Angles of an Electrostatic Gyro 90 


9 


LIST OF 

Table 

Noc 

3.1 

4.1 

6o 1 

7.1a 

7.1b 

7.1c 


TABLES 


Data of the G 10 Two-Degree-of-Freedom Gyroscope 


Coriolis and Centrifugal Acceleration for Three 
Types of Air craft Flying Horizontally at 45° 
Latitude 


Percentage Change of Coning Frequency 




of the 


Moving Gyropendulum for Constant Angular Momentum 
with Respect to Frequency Q ^ ~ cOq of the Stationary 
Gyropendulum 


Error Terms, Affecting the Operation of the Gyropen- 
dulum and Causing Each a Computational Ground Speed 
Error of 1 km/h 


Error Terms, Affecting the Gyropendulum' s Rotor Speed 
Control and Causing Each a Computational Ground Speed 
Error of 1 km/h 

Error Terms, Affecting the Operation of the Azimuth 
Gyro and Causing Each a Computational Ground Speed 
Error of 1 km/h 


1 . NOTATION 


m/sec' 




m/sec‘ 


m/sec ‘ 


m/sec‘ 


c , c 
x’ y 


dyne - cm 
rad/sec 


f(t-t) 


m/sec‘ 


g 

h 


m/sec 
m or km 


Acceleration of the gyro vd.th respect to 
an earth fixed point 

acceleration of the gyro vd.th respect to 
inertial space 

Coriolis and centrifugal acceleration due 
to horizontal flight over the rotating 
earth (s. eq, (A21a) to (A25a)) 

Coriolis acceleration due to vertical 
velocity over the rotating earth 
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damping coefficient of the x, y gimbal 
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nongravitational specific force which the 
gyro exerts on its support (s. eq. (A19) 
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delayed step function (s. eq. (4.16)) 
gravity 

altitude of the vehicle above the earth's 
surface 


1 

m 



r 

s 

t 

u 


rad or deg 
gm 

cm, mm 

l/sec 

sec 


longitude angle 

mass of the pendulosity 

operator 

lever arm of pendulosity 

Laplace operator 

time 

abbreviation (s. eq. (B15) and (Bl6)) 
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V 


C 


F , F 
x’ y 


dyne - cm 
Volt 


F' 


Volt 

Volt 


G 

H 


m/sec 

c 

Rin — cin 
sec 


H' 


gm - cm 
sec 


H' 


gm - cm 
sec 


I , I 

x’ y 


gm - cm 


K = 

L 

M 

P 

R 


|/iR 


rad 

m/sec 

rad or deg 
djme - cm 
km 

m or km 


S 

T .see 

T® = 8^,4 min 
o 


abbreviation (s, eq. (BI5) and (B17)) 
coordinate transformation matrix, e.g. C 

ng 

(eq. (a 8)) for transformation from the 
gyro to the navigational frame 

transfer functions in the gimbal follow-up 
control (s. Fig, 6.2) 

abbrevation for a lead-integral transfer 
function (s, eq. (6.17)) 

gravitational field intensity 

angular momentum of the gyro rotor; in the 
gyropendulum this is supposed to match the 
tuning condition (s. eq. (4.10)) 

angular momenttim of gyropendulum for mat- 
ching the Schuler tuning (s. eq. (3.7a)) 

total angular momentum of the gyropendulum 
(eq. (Al)) 

moment of inertia of the x and y gimbal 
(s. Fig. 6.2) 

gain of the gyropendulum (s. eq, (5 .Ha,b)) 

latitude angle 

torque applied to the gyro 

position 

radius of the trajectory with respect to the 
center of the eeurth 

signal generator (s. Fig. 6.2) 

lead time, integration time (s. eq. (6.17)) 

Schuler period on the surface of the earth 
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V 

m/sec or km/h 

velocity of the vehicle 

(I)AVjj’ 

m/sec 

abbrevations for computational velocity 

(DAVg* 

m/sec 

error terms (s. eq. (C33) to (C35) and 

(i)Av^®' 

m/sec 

(7.13)) 

(i)Av^e' 

m/sec 


6 ( t-T ) 

l/sec 

impulse function delayed vdthT 


rad or deg 

pitch angle, positive, when tall goes 



down 

X 

rad or deg 

celestial longitude (s. eq. (AI 5 )) 


rad or deg 

roll angle, positive, when right wing 
goes down 

4 > 

rad or deg 

yaw angle, positive, when rotation about the 
downward pointing axis is positive; zero for 
northerly moving vehicle 


rad/sec or deg/h 

gyro drift 

= 7.29 

• 10 rad/sec 

earth angular rate 


rad/sec 

auagulau: rate of the navigational frame with 
respect to inertial space 

uj® = Vg/R 

rad/sec 

Schxiler frequency 

w® = 1.235 

“3 

•10 rad/sec 

Schuler frequency on the surface of the earth 

( 

0 

rad, min, sec 

1 

attitude of the spin vector with respect to 
a coordinate freune, earmarked in the 
subscript 

For the azimuth gyro: 


CD 

M 

«• 

<D 

rad, min, sec 

attitude angles of the spin vector with 
respect to the equatorial tangent frame; 
zero when spin axis is pointing northward, 
psurallel to eeirth axis 
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0 + 0 + 
Eo ’ Do 

rad, min, 

sec 

initial misalignment angles 

0y.'e.; ; 

rad , min , 

sec 

attitude angles of the spin vector, measured 
by gimbal-mounted signal generators 


rad 


attitude angles of the spin vector, measured 




by the pickoffs of an electrostatically 
supported gyro (ESG) 

(I) © 

„ 



min,. 

■ • 

inaccuracy of the signal generators 

For the gyropendulum: 



®N’ ®E 

rad, min. 

sTc 

attitude angles of the spin vector with 
respect to the navigational frame; zero,' 
when the spin vector is pointing parallel 
to the vertical, upward for sensor no. I, 
downward for sensor no. II 

0 + 0 + 
No’ Eo 

rad, min. 

sec 

initial value for 0^^ and 0^. 

®No’ ®Eo 

rad, deg. 

sec 

angles for the initial static equilibrixun of 
the gyropendulums (s, eq. (4,6), (4.7) and- 
(4.9)) 

®N1’ ®E1 

rad, deg. 

sec 

angles for the quasi static equilibrium of 

©N2’ ©B2 



the gyropendulums (s. eq. (4.4) and (4,5) or 

rad, deg. 

sec 

(4.24) and (4.25)) 

© 



multiplication node (s. Fig. 4.1) 

T 

sec 


time delay 

Q 

l/sec 


natural frequency of the gyropendulum in- 

2 

eluding the vertical acceleration p h 
(s. eq. (A27)) 


l/sec 


natural frequency of the gyropendulum exclu- 

2 

: ding the. vertical, acceleration p h (s, eq. 




(A30)) : , 

ro 

o o 

l/sec 


initial natural frequency of the gyropendulum 


for matched initial tvining (H^) 


n'lHo 


Subscripts 

a 

b 

g 

h . 

n 

t 

V 

y» z 

E 

I 

N 

D 

0 

1 

I 

II 

Superscripts 

a 

b 


natural frequency of the gyropendulum on the 
moving vehicle for constant angular momentum 
H 

o 

free azimuth frame 
body frame 
gyro frame 
horizontal 

navigational frame 

! • 

equatorial tangent frame 

I 

vertical 

gimbal axes (s. Fig. 3.2«and 3.5) 
east 

inertial freuae 
north 

down, orthogonal to east emd earth axis 
initial value, stationary on the ground 
quasi static equilibrium (s, eq. (4.4) and 

(4.5)) 

gyropendulum no. I, spin vector pointing up 
gyropendulum no. II, spin vector pointing down - 

free azimuth fr^uae 

4 . 

body frame 
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C ' ■ 


cm 

cmd 


coning frequency, natural frequency of the 
gyropendultim disregarding the vertical accel- 
eration p^h 

controlled member, gimbal 
comraeind (torque) 


d 


disturbance (torque), drift 


e 

g 

i 

n 

qs 

s 

s 

T 

* 

Symbols 


earth 

gyro frame 
inertial space 
navigationaj, frame 
quasi static. 

Schuler (frequency, period) 
spin, in connection with 
transpose of a matrix 
computed or measured values 

inaccuracy (preceding the specific term) 


(U) 


1 

6 


A 


uncertainty (preceding the letter M, i.e, 
uncertainty torque) 

unit vector of the angular momentum 

difference ;of corresponding readouts of gyro-; 
pendulum no. I and II, s. eq. (8.16) 

difference from stationary values 

variable in the Laplace domain 

vector 


mean value 


2. INTRODUCTION 

Most of the inertial platforms for navigation with respect to the surface 
of the earth incorporate two or three gyros (two, where two-degree-of-free- 
dom gyros are used; three, where single-degree-of-freedom gyros are used) 
and two accelerometers as inertial sensors. These platform mounted sensors 
are interconnected so that the performance of the platform in a dynamic 
environment is similar to a pendulum with its mass located in the center of 
the earth and with its lever arm equal to the earth radius. This mathematical 
pendulum, if it would be practicable would always indicate the vertical. 
Schuler [1 ] has shown that a physical pendulum or a gyropendulum with the 
same period of T = 8^ min as this mathematical pendulum will indicate the 
vertical if initially aligned with it, even on a moving base. Wrigley [2] 
introduces for the tuning of an instrument to 84 min the expression "Schuler 
Tuning" . 

The history and development of Schuler tuned platforms using gyros and 
accelerometers are well covered in the literatiire l3»4,5,6,8]. A different 
approach to a vertical indicating platform system was carried out by Xstrom 
and Hector [9]i who use instead of three SBF gyros and two acqelerometers, 
two pendulous SDF gyros (implemented by attaching two SDF gyros to two phy- 
sical pendulxims) and one nonpendulous SDF azimuth gyro mounted on a common 
stable platform. As the title [9] indicates, the Schuler pendulum is mate- 
rialized by raising artificially the moment of inertia of the physical pen- 
duliim through the application of a command torque proportional to the angular 
acceleration of the pendulum. The signal proportional to the angular accele- 
ration is generated by the pendulxim-mounted SDF gyro with integral feedback. 
The feedback signal applied to the gyroscope is proportional to the velocity 
of the instrument with respect to inertial space. This gyropendulum, if 
Schuler-tuned will indicate the vertical about the input axes. Koenke [10] 
presents an analysis of this system. 

The system described in this study incorporates also three gyros_copic sen-^ 
sors, two gyropendulums and one azimuth gyro and as such bears resemblance to 
the one developed by Xstrom and Hector. In contrast to the latter. each of the 
three two-degree-of-freedom gyroscopic sensors is base-motion isolated ;ab out. 
its two input axes by means of two gimbals (one gimbal being the gyro case) 
and a follow-up control. The three gyro-gimbal units are mounted on a common 
rigid base which is strapped down to the vehicle. The proper tuning of the two 
gyropendulums when stationary is achieved by a mere center of support and 
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center of gravity separation which is for a modern gyro (e.g. Autonetics GIO) 
in the order of 0.135 nun* One of the two gyropendulums has its spin vector 
pointing upward, the other one downward. 

Similar to a vertical indicating system where the vertical is physically 
indicated by the platform, in the system described herein the vertical can be 
visualized as the mean angle between the spin vectors of the two gyropendu- 
lums. As distinguished from other vertical indicating systems, where velo- 
city is generated by an electric integration process, the two gyropendulums 
indicate physically their velocity with respect to inertial space by the tilt 
angle of their spin vector from the vertical. 

The basic performance of two gyropendulums with their spin vectors pointing 
up and down was first described by Schuler in 1923 [l] who derived the timing 
condition of the sensors for a non-rotating earth. It is shown in this ana- 
lysis that for a vehicle moving over the rotating earth the tuning condition 
is affected by the vertical Coriolis and centrifugal acceleration, the verti- 
cal component of earth rate and the rate of change of latitude. Thus, the 
tuning condition becomes a function of ground speed and position and has to 
be matched by a proper control of the angular momentum of both gyropendulums. 

Wrigley [2], in an interpretation of Schuler's paper, mentions the pro- 
perties of the two sensors as derived by Schuler. Reisch [ 3 ]» iu the descrip- 
tion of the development of his first vertical indicating platform during 
World War II in Germany, mentions Schuler's idea with the two gyropendulums, 
but rejects their possible use as sensors in an inertial navigation system. 

In the discussion of Reisch 's paper, Fischel points out that they might be 
used in such a system. In his paper [8] Fischel describes Schuler's idea, but 
points at the disadvantage that the two gyropendulums indicate only their 
velocity with respect to inertial space. Indeed, to compute herefrom ground 
speed, which is of main interest for terrestrial navigation, earth rate hats 
to be known in magnitude and direction. Schuler proposes to measure this 
vector externally by means of star sighting. In a self contained navigation 
system the direction of earth rate can be stored physically in a gyro, am 
idea which Stratton [l5] also briefly draws attention to, when he discusses 
the behavior of the Schuler tuned gyropenduliun in the spherical gravitationaj. 
field relative to that of the physical pendulum. 

In the navigation system described in this paper the spin vector of the 
azimuth gyro which stores this direction of earth rate is pointing northerly, 
parallel to the earth axis. The magnitude of earth rate is updated in a com- i 
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puter. 

If the azimuth gyro is gimballed as mentioned above, gimbal lock may occur 
when latitude plus pitch or roll angle become 90°, i.e., when the vehicle is 
manoeuvering at high latitudes. This limitation may be circumvented, if 
necessary, by putting additional gimbals between the gyro and the vehicle or 
by replacing the gimballed gyro with a non-gimballed electrostatically 
supported vacuum gyro with direction cosine pattern readout. 

With its sensors mounted on a common base strapped down to the vehicle the- 
described system has advantages of a strapdown system (easy mounting and 
maintainability of the sensors). As. against a strapdown system it. avoids 
error sources due to the dynaunic input to the sensors and the computational 
process. Insofar it has advantages of a platform system, too and can be 
thought of as being a hybrid system. 



3. BASIC PERFORMANCE AND IMPLEMENTATION OF THE NAVIGATION SYSTEM 

In this chapter we will review Schuler's results [l ] with regard to the 
quasi static performamce of two gyropendulums with their spin vectors poin- 
ting up or down, respectively. Furthermore, we will see how an inertial navi- 
gation system can be implemented, using these gyropendulums and an azimuth 
gyro as sensors. 

Schuler [l] considered only a vehicle moving parallel to the surface of the 
earth and neglected side effects resulting from Coriolis accelerations and 
the gyro precession due to the vertical component of earth rate and rate of 
change of longitude. He derived expressions for the quasi static equilibrium 
of the gyros when the torque caused by gyro precession due to the horizontal 
component of earth rate (O)^® cos L) and rate of change of longitude (^cos L) 
and latitude (L) is compensated by a pendulous torque mgr 0 (s, [l], eq. (-A7), 
mr = pendulosity of the gyro, 9= inclination to the vertical, in [l] this is 
called p). Let us introduce the northerly and easterly velocity instead of L 
and I cos L 


Vj^ = R pL 

Vg = R cos L'pl, 

where 


R = R + h 
0 


(3.1) 

(3.2) 

(3.3) 


is the radius of the trajectory with respect to the center of the earth and 


P =■ 


dt 


( 3 .^) 


is an operator. From the equilibrium of the two torques we might obtain the 
following expression for the angles about the north and east axes of the 
gyropendulum for a vehicle which moves with constant velocity parallel to the 
surface of the earth 


^qs _ - H' 
'='nI,II “ mgr 


'^EI,II ^ mgr 


N 


H. 


^ (Rcj-® cos L + Vg), 


(3.5) 

(3.6) 


where the subscripts I, II and the signs + belong to a gyropendulum with its 
spin vector pointing up or down, respectively. 

Schuler has shown in [l], eq. (AlO) that these equations hold even for a 
horizontally accelerated vehicle under the above mentioned assumptions if the 
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gyropenduliun is tuned to 


BSL - i/rz 

H» ■ VR 


Oi 


which is on the surface of the earth 
0 )® = 1.255 * 10 "^ sec“^ 


(3.7a) 


(5.7b) 


or 


T 


= 84.4 min. 


( 3 . 8 ) 


Then we may rewrite eq. (3.5) and (3.6) as 




le 

U) 


cos L + Vg) = + K 


where V.„ is the easterly velocity with respect to inertial space, 
ih 

K = 1/ l/gR 

and on the surface of the earth 

K = 1/ V^R = 1.265 • 10"^ . - 7 — ' = 7.24 ~ 7 ~- . . 
o ' '®o o m/sec km/h 


(3.9) 

( 3 . 10 ) 

( 3 . 11 a) 

( 5 . 11 b) 


The eqs. (5.9) and ( 5 . 10) show us that, the deviation of each spin axis from 

the vertical takes place in the plane normal to the motion of the vehicle 

with respect to inertial space. According to this simplified analysis 

(Chapter 4 represents a more detailed analysis) we obtain for the deviation 

of a gyropendulum which is stationary at the equator (V„ = 0, V.„ = R = 

rJ irj o 

1 , 675 . km/h) 

(3.12a) 

= + 3.37°. ' ( 3 . 12 b) 


For the implementation of this gyropendulum a sensor may be used similar to 
the Autonetics GIQ two-degree-of-freedom gyro. It is described in [ 11 ]; its 
sketch and some numbers are shown in Fig. 3.1. The rotor is supported by a 
hydrodynamic spherical gas bearing providing freedom about all three axes. 
From the gyro data mentioned in Fig. 3.1 and eq. (3.7) it cam be derived that 
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Table 3«1 


DATA OF THE G10 TWO-DSGREE-OF- FREEDOM GYROSCOPE 


(Taken from [ll]) 


Dimensions 


Height : 
Diameter: 
Weight : 


8.^ cm 

8.9 cm 

907 . gm 


Rotor Data 


Moment of Inertia: 


2 

gm-cm 


Mass : 

135 . 

gm 


Angular Rate: 

60 . 

5 240. g 

rps 2 

Angular Momenttim: 

3.6 X 

10 -^ 1.44 X 10 

gm-cm /s 

Drift Data 

Compensated Bias rms: 

0.1 

0.1 

deg/h 

Random Drift rms: 

0.005 

0.02 

deg/h 

Bias Stability rms: 

0.3 

0.3 

deg/h 

Temperature Sensitivity: 

0.005 

0.005 

deg/h/F 

Drift Due to 

Rotor Mass Unbalance: 

0.5 

0.5 ■ 

deg/h/g 

Rotor Structural Compliance: 

0.4 

0.4 

deg/h/g 

Incompressible Effect of _ 

Gas Bearing: 

9.6 

0.4 

deg/h/g 

Compressible Effect of _ 

Gas Bearing: 

1.6 

0.4 

deg/h/g 



Fig. 3»1 A Gas-Film Supported Free-Rotor Gyro 
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Gyro xy- Induction Gyro Case 



Fig, 3*2 Gimballed Gyropendulum 









this gyro ‘is Schuler- tuned if the center of gravity is below the center of 
support by 

r = — ^ = 0.135 mm. (3.13) 

m l/gR 

We will see in the next chapter that the tuning condition is slightly 
different for a gyro on the moving earth. 

The G10 gyro is a null sensor and it must be isolated from the angular mo- 
tion of the vehicle about the two input axes. This might be accomplished by 
putting a gimbal between the gyro case and the vehicle which can rotate with 
respect to the vehicle about the longitudinal axis. The gyro case may rotate 
with respect to the gimbal about the pitch axis. A sketch of this setup is 
shown in Fig. 5*2. We see that the gyro case (shown as innermost gimbal) is 
indeed base-motion isolated about the longitudinal (x-) axis and pitch (y-) 
axis, but is not isolated from the motion of the vehicle about the yaw (z-) 
axis. Since this last mentioned angular motion occurs about the gyro's spin 
axis, it has a negligible effect on the motion of the rotor, especially if 
the coupling torque due to changes of the relative angular rate between the 
case and the rotor (i.e. the viscous aerodynamic torque and the motor torque) 
is small. The relative angular motion between the gyro case and the rotor 
does affect the readout of the gyro's attitude and the command torques, as 
will be seen later. 

On tlie 6Lxis between the vehicle ard the gimbal as well as between the gim- 
bal and the gyro case, two signal generators and two torque generators are 
mounted. The latter ones are connected to their respective gyro pickoffs via 
an electronic network in order to null the pickoff signal between the case 
and the rotor. 

The two si gnal generators read the attitude of the gyropendulum in the 
body frame. The angles 0^ about the longitudinal axis (positive when, the 
vehicle moves positively with respect to the gimbeuL) and 0 about the pitch 

y 

axis (positive, when the gimbal moves positively with respect to the gyro 
case) are related to the angles 0^^ of the gyropendulum about the north axis 
and 0_ about the east axis in the following way (s. eqs. (A45) and ikkS) in 

Hi 

Appendix A) 

0 ^ = vp - 0 ^ cos tp. - 0 g sin (p ( 3 . 1 ^) 

0y = ■& “ ©2 cos tp + ©jj sin <p , ( 3 . 15 ) 
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where we have assumed that the sequence of rotation between the navigational 
fraune (x pointing north, y pointing east, z pointing down) and the body 
frame was first about the down axis through the angle ( for northerly flying 
vehicle (|j = O) , secondly about the pitch axis through the angle •& and 
thirdly about the longitudinal axis through ip , with and •& being small 
angles. 

The eqs. (3.1^) and (3*15) prove Schuler's results: under the assumptions 
mentioned above and using eqs. (3*9) and (3»'10) the sum of the quasi static 
readout of both gyropendulxims indicate the vertical 

‘®xl ‘ 5 . 16 ) 

‘®yl ‘5.17) 

and the difference is proportional to the velocity with respect to space 

^®xl - ^ K 1 ^5.18) 

=6 0yi^ii/2 = K [v.g cos(jJ - sintpj. (3.19) 

Solving for Vj^ and 

\ = -k ‘ «®xl?k (3.20) 

^ ^ ^ ^ (3.21) 

d.6 

where for a known latitude L, the velocity R(u cos L of a point on the siir- 
face of the earth is a known quantity and can be subtracted. The position of 
the vehicle, i.e. latitude and longitude are easily obtainable from V and V„ 

JN £1 

by an integration (s. eqs. (3.l) and (3*2)). In these equations the azimuth 
angle ip is still unknown. The information on the north-direction must be 
obtained externailly based on star sighting, as Schuler proposes it in [ 1 ], 
or it must be stored in a free gyro. In the latter case the resulting navi- 
gational system would be self-contauLned. 

For such a self-contained navigational system another, but non-pendulous , 
GlO-gyro may be used for the azimuth reference. The spin vector of this gyro 
is pointing northerly, parallel to the earth axis. A gimbal between the gyro 
case and the vehicle provides base motion isolation of the gyro about its 
input axea^ as- is shown in Fig. 3. 3. From the two output signals of the signal 
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Gimbal z-Torquer 



Fig;, 3,3 Gimballdd Azimuth Gyro 

(Signal paths from gyro pickoffs to gimbal torquers 
not shown here) 
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generators which are the angles 0 (positive, when the gimbal rotates posi- 

y 

tively with respect to the gyro case) and 0 (positive, when the vehicle 

2S 

rotates positively with respect to the gimbal) one might extract (s. Appen- 
dix B, eqs. (Bl8) and (B19)) latitude and azimuth 


0 a- 

y 


L + ■& cos ij) + sin 4) - 0-j 


( 3 . 22 ) 


0^ a 4^ - ( •& sin 4^ - 'P cos 4^ ) tan L + O^y^cos L (3.23) 

where 0„ and 0^, are error angles due to gyro drift and misalignment, 

Xi D 

Fig. 3.^ shows the signal flow diagrajn for this computation. 

Since latitude is also a result of the integration of the northerly ground 
speed obtained from the two gyropendulums (s. eq. (3.20)), one has in eq. 
(3.22) a redundant information. 

Unfortunately, gimbal lock occurs in this kind of azimuth gyro for 
0^ = - 90°, when the readout of the z signal generator deteriorates (s. eq. 
(BI3)). At 60° latitude this happens, for instance, when the airplane, making 
a turn with 4> = -30°, is heading towards east. In Appendix B, Section B3 it 
is shown that this does not occur, for instance, if an electrostatic gyro 
(described in [l2]) with direction cosine pattern readout is used instead of 
the gimballed gyro. The readout equations of this kind of gyro are (s. eqs. 
(B21) and (B22)) 

01 90®+ L - -S’ cos 4^ - 4>sin 4^ + 0„ (3.24) 

y 

0^ = arc cos [ - sin 4^ cos L - ( vp + 0^ sin4^ ) sin L + 0^^ cos 4> ] . (3.25) 


In the following we will assume that a gimballed gyro is used for azimuth 
reference . 

Compared to an azimuth gyro with its spin vector pointing north but parallel 
to the surface of the earth, in which case 0 provides heading information 
only, eq. (3.23) is more complicated. But with such a gyro a command torque 
has to-be applied- to the gyro for compensating the vertical component of earth 
rate. This causes difficulties insofar as with our proposal the gyro case is 
not base-motion isolated about the spin axis and the output signal of the 
gyro yz-torquer would have to change with the motion of the vehicle. With the 
chosen alignment no command torques have to be applied to the gyro, since the 
gyro spin axis is supposed to stay inertially fixed. 
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From Azimuth 



( Error ) 


yig* SignaJ. Flow Diagraa for the Computation of Azimuth (tj)) and 

Latitude ( LJ j Based on the Readout of the Azimuth Gyro and 
Attitude, Dez^ved from the Gyropenduluma 
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= Rpl 



RU)'* 

Fig. 3.5 Signal FI9W Diagram for the Simplified Computation of Attitude (<Pt ■& ) 1 Ground Speed 

(V„, V_) and Position (L, 1 ), Based on the Quasi-Static Readout of 2 Schuler-Tuned 

N E 

Gyropendulums 







Fig. 5.5 shows the signal flow of the navigational computation for the 
vehicle's attitude ( ip , -d , eqs. ( 3 . 16 ) and (5.17)) with respect to the ver- 
tical, the vehicle's north and east velocity (V^,, Vg, eqs. (3»20) and (5.2l)) 
and finally the vehicle's position, namely longitude S and latitude L. The 
azimuth angle 4^ is computed from the azimuth gyro readout according to 
Fig. 5-^. 



k, PERFORMANCE OF THE GYROPENDULUM IK A DYNAMIC ENVIRONMENT UNDER IDEAL 
TUNING AND COMPENSATION CONDITIONS 

In the previous chapter we have seen that one may obtain from the sum and 
difference of corresponding quasi static readouts of two Schuler tuned gyro- 
pendulums with their respective spin vectors pointing up or down emd one 
azimuth gyro, all the information needed for navigation when moving over the 
earth, namely the vehicle's attitude, azimuth, ground speed and position. In 
this chapter we will look more thoroughly into the dynamics of the gyropendu- 
lums. This allows us to draw certain conclusions on the practicability of a 
navigational system based on these gyropendulums as inertiel. sensors. Since 
all the navigational information, for known azimuth (j) , is already included 
in the quasi static output signals of the gyropenduliims , any superimposed 
oscillation deteriorates this information and has to be suppressed by means 
of control loops and proper initial alignment. 

Eq. (A26) of Appendix A describes the dynamics of a gyropendulum ; Fig. 4.1 

is the corresponding signal flow diagram. It shows that any torque about the 

emd 

east axis caused by the torquer (M^ ) or a northerly acceleration fjj causes 

the gyropendulum to tilt about the north axis with respect to inertial space 
Since the earth is also moving with respect to inertial space and 
since the motion of the gyropendulum is only affected by the tilt from the 
vertical, the angle enters the loop at the node 1N. The tilt from the 

vertical 0^^ = causes the gravity and the vertical component of 

earth rate to generate a torque about the north axis at the node gN, which 
stimulates the motion about the east axis. 

From the left hand side of eq. (A26) we see that in general the spin vec- 
tor of the gyropendulum describes an undeunped motion - actually a coning 
motion as we will see later - the frequency (s. eq. (A27)) of which 

SS (1 + _ ^) 1 + -^) ein l (4.1) 

depends on the vertical acceleration of the vehicle and the vertical compo- 
nent of earth rate and rate of change of longitude or easterly velocity, 
respectively. For the same parameters m, r and H of the gyropendulum, the 
frequency of the instrument with upward pointing spin vector is slightly 
different from the one with downward pointing spin vector. 

2 

The vertical acceleration p h han zero mean indeed for a horizonteLLly 
flying vehicle, but short time peaks due to verticail gusts might very well 
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Pendulous Gyrowheel (North and East Axis) as 
Torque Summing Member 

Relative Motion (North and East Axis) betveen 
Gyropendulum and Navigational Freuae 

Multiplication Nodes 

Mechanical Signals Magnetic Signals 




4«1 Signal Flow Diagram for a Two*Degree*of«Freedom 
Gyropendulum with Vertical Spin Axis 











exceed 1g. Since these peaks last much less than the period of the gyropend- 
ulum, which is in the order of 8^ min, we have, in order to estimate its 
effects upon the motion, approximated the acceleration of one gust by am im- 
pulse function occuring at the time T (s. eq. (A28)) 

p^h aph*5(t-T). (4.2) 

Since this form differs from zero only for a short time interval at t = T , 
it will be multiplied in eq. (A26) with the constant angles 0^^ (t ) or 0^ (x), 
respectively. This shows us that the short time vertical acceleration causes 
an additional forcing function which enters the right hamd side of the equa- 
tion (s. eq. (A56)) while leaving the frequency Q on the left hand side 
(s. eq. (A30)) given by 


n® = 


H 


(1 - i ((0 


le 


g 


Rc osL.' 


) sin L 


(4.3) 


Based on this slowly varying frequency on the left hamd side ofeq. (A26), 
we have extracted the quasi static equilibrium, i.e. , the equilibrium between 
the initial torques on the right hamd side (subscript ”0**) and the slowly 
varying "spring torques" Hi) 0^^ on the left hand side. This quasi static 
equilibrium is (s. eqs. (A34) and (A35)) 




(4.4) 



M^^/H cos L 


(4.5) 


These values are essential, for the initial alignment of the gyropendulum. 
Assiaming that a commamd torque is applied to it for deunping out the oscilla- 
tion during initial self aLlignment mode, it will always settle to this equi- 
librium which at 45° latitude is for = 0 amd ~ UJ® 



(4.6) 


0 gQ ~ + 2.39° • 


(4.7) 


If the disturbance torque about the north aixis, expressed, like the gyro 
drift, in degrees per hour, has the following magnitude 
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(4.8) 



1 MERU = 0.015 deg/hr 


the gyropendulvun vd.ll settle with an offset angle about the north axis from 
the vertical of 

0 ,, a 12.2 sec . (4.9) 

No 


As was pointed out in the previous chapter, the quasi static output signal 
of the gyropendulum mounted oh a moving vehicle conted.ns all necessary navi- 
gational information for knovm azimuth. Any superposed oscillation deterio- 
rates this information. The remaining part of this chapter deals vd.th the 
requirements to suppress the oscillatory motion by means of control loops and 
proper initial alignment of the spin axis. 

From eq. (A37) to (A39) in Appendix A the follovd.ng may be concluded: 
tuning the gyropendulum to 


_H 

mr R 



(4.10) 


causes in eq. (A 38 ) and (A39) the parenthesis of the first term on the right 
hand side to cancel vd.th the determinant A and the second. term to become 
zero. Physically this mesuis that the gyropendulum vd.ll indicate the north and 
east velocity and vd.thout time lag and oscillation and that it becomes 
insensitive to the horizontal accelerations pAV = sAV. Since we have not spe- 
cified in Appendix A the transient of AV^^ and except for continuity 

(AV„^ =AV.^ = 0 ), this result holds for auiy continuous motion of the 

No xEo 

vehicle. 

The tuning condition (4.10) can be easily derived from the signal flow dia- 
gram in Fig. 4.1, keeping in mind that, caused by the earth's geometry, amy 
northerly acceleration, for instance, is coupled vd.th an amgulao* rotation 
about the east axis ( ^ jj = = RsWg^). Proper tuning reqtiires that the 

signal flow, originating in for instance, is not circulating in the closed 
loop but ends in the point IE. This means that the northerly accelerated gyro- 
pendulum tilts about the north axis (0jj = mr f^/sE = mr/H) but tracks the 
vertical about the east axis (0„ = 0 or 0i® =0^^) . 

ill a* 

Before we continue in the derivation of the performance equations of the 
gyropendulum, let us first have a look at the Coriolis and centrifugaG. accel- 
eration for three types of aircraft: a private plane or helicopter, a sub- 




34 



sonic jet and a supersonic transport, flying at different speeds in norther- 
ly, north-easterly and easterly directions. For these aircraft, equations 

(A21a) to (A23a) have been evaluated. As we see from Table 4.1 these accelera- 

-2 

tions are very small, the largest value being b^^ = b^ = 1.5 * 10 g for the 
SST flying at 2^500 km/h in an easterly direction. But it can be easily shown 
that already the horizontal Coriolis acceleration of the private plane or 
helicopter has a serious effect on inertial navigation. In order to prove 
this physically, let us look at a normal pendulum subject to a constant hori- 
zontal Coriolis acceleration b^ in the gravitational field. The ratio t)^/g 
is the value of the steady state angle which the pendulum indicates with 
respect to the vertical. The same is true of the gyropendulum and since 
1 s^ tilt angle of the gyropendulum corresponds to 0.l4 km/h ground speed, 

(s. eq. (3«1'l)) the indicated ground speed error caused by the horizontal 
Coriolis acceleration is 

AV = 2.9 * 10^ km/h . (4.11) 


Table 4«1 Coriolis and Centrifugal Acceleration for Three Types of Aircraft 
Flying Horizontally at 45° Latitude 


No. 

Flight 

Direction 

Components of 
Coriolis and 
Centrifugal 
Acceleration 

Aircraft 

Private 
Plane , 
Helicopter 

Subsonic 

Jet 

Supersonic 

Transport 




V = 250 km/h 

V = 900 km/h 

V = 2; 500 km/h 




h = 1 km 

h = 10 km 

h = 20 km 

1 


bjj/g 

0 

0 

0 

2 

North 

bg/g 

5.65 * 10"^ 

1.52 • 10“^ 

5.65 • 10”^ 

3 

(V =V V =0) 

\/s 

7.70 • 10"^ 

1.00 • 10"^ 

7.70 • io"5 

4 


^n/s 

5.50 • 10"^ 

2.55 • io~^ 

9.10 • 10“^ 

5 

North-East 

bp/g 

2.98 * 10"^ 

1.43 • 10“^ 

6.45 • 10"^ 

6 

(Vjj=Vg=v/y^ 

Zi 

\/s 

5.90 • 10"^ 

2.86 • 10“^ 

1.19 • 10"^ 

7 

East 

bjj/ S “ by/ g 

8.15 * 10"^ 

3.62 • 10“^ 

1.50 • 10"^ 

8 

V''> 

bg/g 

0 

0 

0 
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From Table 4,1 we see that for a private plane flying easterly, this 
error amounts to AV = 23 km/h! Hence these error terms have to be very well 
taken care of. 

In order to judge the usability of the gyropendulum as an inertial sensor 
we will now derive from eq. (A38) and (A39) the performance equations under 
the following aesumpticaiB which have to be matched in the real world: 

a) the torques caused by the Coriolis and centrifugal accelerations of a 
horizontally flying vehicle are properly compensated, i.e. 

(4.12) 


,,cmd , 

Ml = -mr b 


„cmd , 


(4.13) 


b) the gyropendulum is properly tuned according to eq. (4.10), 

c) the disturbance torques acting on the gyropendultim do not change from 
their static values, i.e.AM^^ = AM^ = 0. 

With these assumptions we obtaiin for the displacement angles of the gyro- 
pendulum about the north and east axis from their quasi static equilibrium 
in eq. (4,4) and (4,5) in the time domain: 

A6 jt = + cos + A0j^ sin fi^t 


N 


No 


So 


siph f (t -T) -T— "©V (^)Vos(n" (t -X)) 

{ x )^ s ± xi ( q ^ (t - I ))] 


(4.14) 


^ ^\e ^^®Eo cos fi'^t iA0j^ sin fi^'t 


mr u ~ 

- — Ph, 1 ( 


-®E -t))] 


(4.15) 


where 


f (t - X ) = 


0 for X > t 


1 for X < t . 


When the gyropendulma^ mounted on a stationary vehicle (AV = AV.„ ^ 

N 

is initially aligned with the verticaOL, i.e.A0^ = -0„ and A0^ = 

’ No No Eo 


(4.16) 

: ph = 0), 
Eo 
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(s. eq. (4,9) and (4.7))» it oscillates in a coning motion about 0^^^ and 
as shovm in Fig. 4.2,: 



Looking 


North 




ng. 4.2 Coning MoUon of n Oyropoadnlim with Spin Vootmr Pointing Up 
or Downt BoapoetiToly 

Let us now turn our attention to the angles 0^^ (^) and 0^, (t) on the right 
hand side of eq. (4.14) and (4.15). The coefficient of these two angles is 
at the time t > X for ph = 1 m/sec 

|£ ph = a; Sh — » i.3 . 10"^, (4.17) 

® Ruf 

o 

having kept in mind eq. (4,10). This may be neglected compared with 1 on the 
left hamd side of eq. (4,l4) amd (4.15) even for higher vertical velocities. 
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Thus we can neglect 0j^ (x) and 0^, (x) on the right hand side as compared 
with the left hamd side, amd for t = X we obtadLn for a perfectly aligned 
gyropendulum neglecting the gyro drift in eq. ik.k) and (4.5) 5 






(x) 




- /(U^^cosLq 
1? 


nir 

eT 




(x) 



(4.18) 


(4.19) 


Introducing these into the right hand side of eq. (4.l4) and (4.15) we see 

2 

that the vertical acceleration p h, which we have assumed to be am impulse 
function at the time x , does not stimulate am oscillation, provided that 
ground speed remains constamt during ascent or descent. It cam be easily 
shown that oscillations due to changes in groimd speed are negligible 
(s. eq. (C23) of Appendix C). 

For the total angle of the gyropendulum about the north and east aixes we 
may write 


0 : 


N 


= ®N1 - 


F 


— mr ,, 

4. ..M }T 

H N 




v.„ 

^ mr lE 


= 


®E1 ^ ^®E 


- mr ^o 


mr N , 
- ^ — ph 




f (t - X ) + A0^^ 


cos 




t i A0+ 
No 


sxn 


Q 


' t. 
(4.21) 


The resTilt that the vertical acceleration does not stimulate am oscillation 
but yields only am angle displacement proportional to the product ph or 
Vjj ph respectively, is really surprising amd we have to find a physical ex- 
planation for it matking use of the signal flow diagram in Fig. 4.1. Assuming 
that we have a perfect and aeeurately aligned gyropendulum (M*^ = 0, A0 ^ = O) 
and that the vehicle is flying eastward, we have prior to the vertical, accel- 
eration 0JJ = 0, 0J. = + ^ \e* vertical acceleration p^h = s^h occurs 

there will be an additionaJ. torque with the magnitude -mr ©tt at the torque 
summing member marked gE in Fig. 4.1. This torque is integrated in the gyro- 
pendulum by i 1/*H, resulting in an offset angle A0^j about the north axis, 
c 

The torque -HQ A0jv catxsed by this offset angle arrives at the torque summing 


k. 
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member gN in Fig, 4.1. After the vertical acceleration, the vehicle has a 
vertical velocity sh which causes an easterly Coriolis acceleration sh V. „/R 
resulting also in, a torque at the torque summing member gN. If both torques 
cancel each other, the signal generated by the vertical acceleration will 
not circulate in the loop and hence will not stimulate an oscillation. The 
mathematical condition for this is ^ 


(- mr s^h 0g) (i 1/Hs) (- HQ^) = - mr sh or 




^iE coincides with the assumption stated above. 


After this consideration of the signal flow the offset angle 4as alsfo 

become understandable. It is the result of the integrated torque caused by 
the vehicle acceleration and we find 


A0- (ph) ^ i mr 


N 


sH 


0 - mr 

E ■ H 


!ii 

rq'^ 


ph 


(4.22) 


which proves the results in eq. (4.20). 

Since the velpcity is computed from the difference of corresponding read- 
outs of both gyropendulums with their spin vectors pointing up and down, this 
offset due to vertical velocity does not affect directly the velocity compu- 
tation. It does affect the attitude computation which is derived from the sum 
of corresponding readouts and via the circumference of the azimuth compu- 
tation it does also affect the velocity computation (s. eqs. (5«23),, (3.»20) ^ 
and (3.21)). But its effect on attitude is not too weighty. For an easterly .. 
flying SST (V^. = 2,500 km/h, = 3,685 km/h at L =. 45°), the offset angle 
about the north or pitch axis respectively for ph = 1 m/sec and with 
l/(RflP) = mr/H ~ K (s. eqs. (4.10) and (4.26)) is 


ph ■ ph iE 


3.4 


sec 

m/sec 


(4.23) 


It will be shown in Chapter 7»2 that its effect on the velocity computation 
is not too significant either. F\irthermore , to achieve higher accuracy, this 
offset angle can be easily biased out by defining as the new quasi static 
equilibrium 
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®N2 = ® 


N1 


A0. 


N 


*Slo mr \e 
+ T 


HQ" 


H 


RQ" 


(4.24) 


®E2 = ®E1 -"^®E 

_ Mg^H ; 0)^^ COB L 

Q° 


mr 

® RQ® 


ph 


(4.25) 


(0jj^, 0g>|» s. eq. (4.4), (4.5)). In an airplane, for instance, a signal for 
the vertical velocity ph is readily available and the compensation can be 
easily accomplished. 

In the follovdng mathematical derivation we will neglect the effect of 
vertical velocity on attitude and velocity and leave the proof for the vali- 
dity of this assumption to Chapter 7.2. 


Eqs. (4.20) and (4.21) show us that under the previously mentioned asstimp- 
tions a to c (compensation for the Coriolis acceleration in horizontal flight, 
constant disturbance torques, perfect tuning) the gyropendulum when initially 
perfectly settled to its equilibrium (^®^^ = indicates indeed the 

velocity of the vehicle with respect to space as stated by Schuler [l] 

(s. eqs. (5«9) and (3»10)), with only the gains seeming to be different. In 
order to examine this, we solve eqs. (4.3) and (4.10) for mr/H, introduce 
this into eqs. (4.20) and (4.21) and calculate the difference and the sum of 
corresponding angles of the two gyropendulums . As we recollect, this is nee- 
ded for velocity smd attitude computation. Keeping in mind that the vertical 
Coriolis acceleration is a small quantity (s. Table 4,1) we find from eqs. 
(4.3), (4.10) and (3.11) 


mr 

H 


» K [d + I ^ * 


RcosL 


) sin L ] 


(4.26) 


where we have neglected in the square root ofK (s. eq. (3*11)) 
r ie ^E 

Kh) + g^^^) (sin/2g) J compared with l/gR. For « 0 and L = 45® 
this term is 1o^ times smaller* 

The angles of the gyropendulum ' s spin vector about the north and east axes 
from the quasi static equilibrium in eq, (4,4) and (4.5) are thus 
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A0„ =6„ - e 


N 




V 

_1_ (uji® + — J ) sin L ] AV-. 
^ RcosL J ^ 


AQ+ 


A0. 


®E ■ ®E1 


= IkI (1 t 7 (W 


le 


2o) 


He osL 


) sin L ] AV_.^ '*' ^®Eo ^ 

(4.28) 


+ A0^ sin n ''t . 
No 


Similar to these equations, it is convenient for the following evaluations to 

have the quasi- static equilibrium of eqs. (4.4) and (4.5) biased out from the 

readout angles 0 and 0 . This may be done according to eqs. (A45) and (a46) 

X y 

in the following way; 

A0X ~ ®x ®N1 ^ ®E1 = »p - A0j^ cos - A0g sin ip (4.29) 

A0y ' ®y ®E1 ” ®N1 ^ = -& - cos (p + A0JJ sin ip . (4.30) 

The navigational computation is based on the difference and sum of these 
readout angles of the two gyropendulums. The computer assumes the unknown 
misalignment angles and from the quasi- static equilibrium to be 

zero and relies on the following basic equations; 

(A0^j - A0^3 -j)/ 2 = K (1 + by2g)(Av* costp + AV^*^ sin 4» ) (4.51) 

(4.32) 

(4.33) 

(4.34) 

where the asterisk denotes a computed value as agsiinst the true value. 

The eqs. (4.31) and (4.32) show us that the gain relating the velocity with 
respect to inertial space and the readout differences is modified from the 
gain in eq. (3.11) by the vertical Coriolis emd centrifugal acceleration and 
we may write: 


(A0yj - A0yjj)/2 = K (1 + by2g)(AV*j. cos(p - AV* sin ip ) 

(A0^j + A0 ^jj)/ 2 = 9*- — tan L (AVj^ cos tp + AV*^ sin ip ) 

^2 

(A0 + A0 ttV2 = ■&*- •=— V.„ tan L (AvT- cosip- AV* sin tp ) , 

yl yll 2 lE lE N ^ ’ 


K* = K (1 + by2g) . 


(4.35) 
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Now we solve eqs. (^.5l) to (^.3^) for the vehicle's velocity and attitude 
and find for the computation carried out in the computer: 


= ''n = - “xll> '=‘>='1' - - ^e „) Sin4,] 

ciK. 


(4.36) 


Avf_ = V* + Rof® (cos L*- cos L ) = 
xE E O 


'p*= 


^*= 



3K’ ‘ 

- A0yij) cos 4) 

+ (AO ^ 
xl 

- AO. 



'*®xl - 

2 

K v.„ 

lE 

tan L 


2 * 

2 

2K' 




(A0 _ - A0 ) 
yi yii 


tan L 


2 ■ 

2 

2K' 



(4.37) 


(4.38) 


(4.39) 


Fig. 4.3 shows the signal flow diagram for the navigational computation of 
the vehicle's attitude (‘P, ), groxmd speed (V„, V„) amd position il. L) . 

The inputs of this diagram are the output signals of the gyropendulums 

II’ ®yl II^’ azimuth angle tj) derived from the output signal of the 
azimuth gyro, the altitude derived from the altimeter for updating R and 
finally the quasi static equilibrium angles tt ^.nd tt of the two 

I X f X X E n X y XX 

gyropendulums. For continuously updated earth rate component amd continuously 
updated disttirbance torques based on the model = f (acceleration, tempe- 
rature, etc) which is derived from laboratory tests, the latter angles are 
known quantities. If the disturbance torques are unknown, they can initially 
be measured when the vertical amd true north are measured externally. We will 
consider the alignment problem in Chapter 8. 
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Pig. 4.3 Signal Plow Diagram for tho Exact Computation of Attitude 
( « '&)« Ground Speed (Pg* V^) xnd Poaltlon (L, 1), Based 

on the Readout of 2 Properly Tuned Qyropendulums 






























5. PERFORMANCE OF THE AZIMUTH GYRO 


Chapter 3 has shown us that a self-contained navigational system may be 
implemented using two TDF gyropendulums and one TDF azimuth gyro. We have 
looked into the ideal performance of the gyropendulums in Chapter h and will 
now examine the performance of the azimuth gyro whose spin vector points 
northerly parallel to the earth axis. 


The basic equations for this gyro in the so called "equatorial tangent 
coordinate frame" (s. Fig. B1 ; subscript "t"; x or N' pointing north, parallel 
to the earth axis, y or E pointing east, z or D pointing down, orthogonal to 
X and y) are derived in Appendix B. 

From the determinant A in eq. (B6) 

A= s^ + ^ (5.1) 

where the frequency is from eq. (B2) 


0 )^^ =0)^® + Vg/(R cos L) (5.2) 

we conclude that this gyro describes an undamped 2k hr mode in a stationary 
vehicle. This is quite obvious if we consider that this non-pendulous gyro 
ideally keeps its attitude with respect to space, i.e. if the spin vector 
does not point exactly northward, it will describe a circular motion aroiind 
true north with a 2h hr period. This oscillation of the displacement angles 
0g and 0JJ about the east and down axis may be derived from eqs. (B4) and 
(B5) as 


0 ^ = 


— 7-r + (0„ rr 1 cos w t 

d 

/q+ ^E \ . , it . 

+ (0„ + — rr ) sxn 0) t 

V Do jjXt j 


(5.3) 


®D = 


‘^E \ it 

+ — rr cos 0 ) t 
Do 


+ (®^o *:h) 

-(®Eo-:it)- 


sin t 


(5.4) 


The errors caused by initial misalignment need no further comment. The ones 
caused by the drift terms oscillate about the static equilibriixm when the 
uncertainty torque Ho) is compensated by the gyroscopic torque 0Ho)^^ which 
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is the component of earth rate experienced by the input axis when the corres 
ponding spin vector is misaligned. It is interesting to note that this error 
angle due to constant gyro drift is greater for a westerly (V_,<0) than for 
an easterly flying vehicle. 

Compared to an azimuth gyro with its spin vector pointing north but 
parallel to the surface of the earth, the chosen orientation seems to be 
more advantageous because the angular rate along the spin axis with respect 
to inertial space, which balances the gyro drift, is the maximum possible 
thus minimizing in eqs. (5*3) and (5*^) the error angles and 

But if we consider the nature of the output signal of the z signal generator 
mounted on the gimbal (s. Fig. 3*3 and eq. (B19)) 

®D 

0 « 4* “ ^ 4^ - cos 4^ ) tan L + ^ ( 5»5) 

we see that the effect of the error angle Oj^ or gyro drift, respectively, on 
the azimuth computation increases with growing latitude like in a conventio- 
nal platform system (s. [13]» eq. (2.45)). 


The output signal 0 of the y signal generators may be used for redundant 

y, ‘ 

measurement of latitude (s. eq. (Bl8)) 

0 - L + ■& cos (|) + vp sin 4^ , - 0- (5»6) 

y ^ 

where and <43 are derived from the two gyropendulums. 
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6. DESIGN CONSIDERATIONS OF THE CONTROL LOOPS FOR THE GYROPENDFLUM AND THE 
AZIMUTH GYRO 

Chapter 4 has shown us that the gyropendulum must be tuned in order to 
indicate the groiind speed without time lag or oscillation (s, eq. (4.10)). 
Since this tuning condition changes with groimd speed and latitude, a loop 
has to be designed for the control of the angular momentum. 

For the gyropendulum and the azimuth gyro, a gimbal servo system has to be 
designed. Finally command torques have to be applied to the gyropendulum in 
order to compensate the torques due to the Coriolis and centrifugal accel- 
erations in the horizontal flight (s, eqs. (4,12) and (4,13))* We will treat 
these control loops in this chapter. 


6.1 Tuning-Control of the Gyropendulum 


The tuning condition (4.10), which is a function of ground speed and lati- 
tude, has to be matched by the proper control of the angular momentiun. In 
order to appreciate the urgency of the timing control, let us first examine 
how strongly the frequency fi® depends on the flight condition if the angular 
momentum is kept constant at that magnitude H^ which provided tuning of the 
gyropendulum on a stationary basis. From eq. (4,26) we find for H^ 


H 

o 



U)^®sin L 
2 o 


where H' is 
o 


H' 

o 

ie j T 
1 I tJL) sin Lg 


2^0 


( 6 . 1 ) 


h' 




®o • 


( 6 . 2 ) 


So according to eq, (3»7a) achieves the Schuler tuning of the gyropendulum 
on the ground which corresponds to the exact tuning condition for nonrotating 
earth. 

With eqs. (6.l) and (4.3) we find for the coning frequencyfi^ on the ground 
^ ~ • (6.3) 

0,02 O 
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Assuming that the actual frequency for constant angular momentum is 


Q 


Ho 




H. 


(6. if) 


furthermore making the following assumptions for R and g 


R = R + h 
o 




g = 


(1 + h/R ) 


2 ’ 


(6.5) 

( 6 . 6 ) 


we find from eq. (if«3) for the change in frequency for constant angular 
momentum 




Ho 




R 


g 






cos L + 
o 


Rcos^L 


) AL + ~ tan L 
R o 


]. 


(6.7) 


where v;e have made use of the approximations sin L k sin + AL cos and 

tan L =: tan L + AL/cos^ L . 

o o 

Table 6,1 shows us that this change in frequency exceeds 11 % for the east- 
erly flying 3ST, 

The change of angular momentum with changing flight condition for matched 
tuning can be easily derived from eq. (4.26) • V^ith 


II = H + Ah 

o 


( 6 . 8 ) 


and l/K ~ (l - AH/H^)/H^ we obtain for this change in angular momentum ex- 
pressed in terms of the frequency change for constant angular momentum 
.01 


Afi 

~ 1 
H ~ 2 

o U) 


H 


o ^ 
s 2 R 


(6.9) 


The second term on the right hand side (h/2R ) takes into account that 


according'to eqs. (4,3) and (4.10) the tuning condition not only depends on 
^ gravity as does the frequency for constant H but also on the radius R of the 
flight path with respect to the center of the earth. 
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Table 6 . 1 Percentage Change of Coning Frequency of the Moving Gyro- 

pendulum for Constant Angular Momentum with Respect to Frequency 
X 00 ® of the Stationary Gyropendulxim 


No. 

Frequency 

Chainge 

Aircraft 

Private 
Plane , 
Helicopter 

Subsonic 

Jet 

Supersonic 

Transport 

V = 250 km/h 
h = 1 km 

V = 900 km/h 
h = 10 km 

V = 2,500 km/h 
h = 20 km 

1 

Afi°(h)/U)® 

o 

- 0.03 

- 0.51 

- 0.63 

2 

AQ^'Cb )/o)® 

V 0 

- 0.08 

- 0.36 

- 1.50 

3 

Ajf(AL)/o)® (V„ = 0) 
o hi 

- 0.57 

- 0.37 

- 0.37 

4 

Afi'^ (AL) Ao® (v„ = V) 

O i!i 

- 0.52 

i 0.92 

- 1.91 

5 

AQ®(V„)/oj® (AL = 0 ) 

th O 

- 0.88 

- 3.15 

- 8.75 

6 

Total: AqVo)® 

- 0.11 i 1.40 

- 0.67 - 3.52 

- 2.13 - 9.12 


According to [11], in the G10 gyro the rotor speed is optically measured 
and can be controlled by means of supply voltage control. The functional dia- 
gram of the gyropendulum in Fig. 6.1 shows this control input. 

At the end of this consideration let discuss the change of the coning 
frequency for matched tuning in terms of the already known frequency change 
for constant angular momentum H^. We write eq. (^.10) in the following form 


=qI + aq'" 


H + AH 
o 

mr R ( 1 +h/R ) 
o o 


( 6 . 10 ) 


from which we may derive the approximation 

AnC 

. Ah ^ ^ I o h 

~”o ■ ">0 ' 30,= ■ ' 

o o . 


Except for the altitude term, the magnitude of this frequency change is half 
the amount of the values listed in Table 6.1. 
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6.2 Gimbal Control and Command Torques 

The principle of gimbal control for the gyropendulum was outlined in 
Chapter 3 ^ind Fig. 3*2. The output signals of the x and y pickoffs in the 
gyro case are fed via an electronic network to the x and y torquers of the 
gimbals . 

The signal flow diagram of the gyropendulum including the gimbal follow-up 
control is shown in Fig. 6.2. The central part of this diagram describes the 
signal flow in the sensor and is the seune as in Fig. 4.1. In order to under- 
stand the kinematics of the different pickoffs namely between the gyro 

(rotor, "g”) and the controlled member (inner gimbal or gyro case "cm") and 
signal generators between the inner and outer gimbal or outer gimbal 

^ » y 

and vehicle ( "b")^ respectively in the gimbal follow-up control, let us first 
have a look at the coordinate frames and the angles between them. (We omit 
in the following the subscripts x and y, since the signal flow is identical 
for both loops). Fig. 6.3 shows one representative axis from each of the 



Fig. 6.3 Angles between Coordinate Frames 

5 frames involved, namely the inertial frame ("i"), the free azimuth frame 
("a", this fraime moves with the vehicle about the vertical z axis with the 
angle ip , but stays horizontal about the x and y eixis), the body frame ("b"), 
the gimbal fraune ("cm" according to "controlled member"), and the gyropend- 
ulum frame ("g"). For the gimbal follow-up control we direct our attention 
first to the angle i.e., the output signal of the gyro pickoff 

This angle shall be driven to zero by the control loop, because from 


gem _ acm ag 

0 - 0 - 0 


( 6 . 12 ) 
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we conclude that then 0 = 0 i.e., the gimbal follows the gyropendulum 

3.C131 

exactly. What are these angles 0 and 0 °? 

The most convenient reference frame for measuring the attitude of the gyro- 
pendulum is the navigational frame, where the readout angles are 0j^ and 0g. 
Since the pickoffs of the gyropendulum move with the vehicle about the ver- 
tical axis, we take the free azimuth frame as a reference frame for descri- 
bing the attitude of the gyropendulum. The angles in this frame are 


0 = 0„ costl) + 0- sin (6.13) 

X N ^ E 

0 = 0„ cos - 0„ sin 4^ (6.l4) 

y k w 


which may be derived from eqs. (A4-5) and (A46) for -9 = ip = 0. Hemember the 
free azimuth frame and the body frame coincide when the vehicle flies hori- 
zontally. 

Let us now have a look at the gimbal attitude 0^°™ with respect to the free 
azimuth fraime. Any torque acting on the gimbal causes it to rotate through 
with respect to inertial, space. Since the free azimuth frame is moving 
with respect to inertial space too, only the difference 

acffl icm ia 

0 =0-0 (6.15) 


effects the readout equation (6.12) of the gyro pickoff 

The gimbal attitude 0^*^® with respect to the free azimuth frame is read out 

by the gimbal signal generator 3*^®^ whose output is the angle But since 

the vehicle moves relative to the free azimuth frame through the attitude 

angles ip and -9 , the readout angle is equal to the difference of the vehicle's 

attitude and the gimbal attitude from the free azimuth frame. We have seen 

eq. (6.12) that for 0®*^"” = 0, the gimbal attitude is equal to the gyro- 

pendulum's attitude 0^® specified in eqs. (6.13) and (6.l4). So in this case, 
cinb 

the angle 0 is truly the angle we seek and which we have derived in eqs. 
(a 45) and (A46). 

The output 0®^™ enters the electronic network and torque generator 
(s. Fig. 3*2 or in Fig. 6.2 the boxes [(l + Ts)/T-b] * F (s)) as an electrical 

signal ( ) and is transformed into a magnetic signal ( ). The transfer 

function for the output of a y- gimbal signal generator say is 
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fe F F' - -&(! + F F' ) 

2 . _ „ I y y y- y 


0 = ■ 

y I s" + c 8 + F F' 


where 


- 0 I + M ], 

y y y ^ 


(6.16) 


F' = 


1 + Ts 
Ts 


(6.17) 


For nearly frictionless gimbal bearings (c-*-0) this simplifies to 

y 


0^ 5S - •& + 


I a + F F' 

y y 


[q p F* - 0^^ I s^ + 

l V V 


M 1. 

y^ 


(6.18) 


For low frequencies (s-*-.0), the output of the gimbal signal generators 
is the desired one s: - -& + especially since according to eq. (6.17) 

F' then behaves like an integrator. For high frequencies F' approaches unity 
and F can be designed to obtain a damped network by putting a proper lead 
into it. 

The performance of the network can be improved if the box F is biased with 

the negative command torque (- M°™*^/F) . This torque is a known quantity 

(s. eqs. (4.12), (4.13), (A49) and (A50)). Then changes in the disturbance 
d 

torque only cause short term errors. 

The design philosophy of gimbal follow-up control is very similar to that 
of the control loop for a servo test table for testing gyroscopes. The latter 
is treated in more detail by LUbeck and the author in [l4]. Appendix 7, and 
it worked satisfactorily, even in the presence of high uncertainty torques. 

We omit the discussion of the follow-up control of the azimuth gyro which 
is very similar to that of the gyropendulum. 
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7. SII'tPLIFIED ERROR ANALYSIS 


The requirements for accurate operation of a self-contained inertial navi- 
gational system using two gyropendulums and one azimuth gyro as sensors have 
been found in Chapters 4 and 5* The following errors in the operation of the 
gyropendulum limit the accuracy of the system: 

a) tuning errors in eq. (4,10), 

b) initial alignment errors (I) and (I) from the equilibrium in eqs, 

(4.4) and (4.5), 

c) errors in the compensation of the biases due to disturbance torques 

(I) and (l) and to horizontal earth rate (l)(U)^® cos 

(s. eqs. (4.4), (4.5), (4.29), (4.30) and Fig. 4.3), 

d) errors in the command torques (I) and (I) for compensating 

Coriolis and centrifugal accelerations due to horizontal velocity 
(s. eqs. (4.12) and (4.13)), 

d d 

e) errors due to the change of the disturbance torques A M^j and AKg and 

f) readout errors (I) 0 and (I) 0 . 

X y 

The system's accuracy is similarly affected by the following errors in the 
operation of the azimuth gyro; 

g) alignment errors s-nd (1)0^^^ , 

h) gyro drift and u)^ and 

i) readout errors (l)0 and (l)0 . 

Y ^ 

For the three types of eiircraft, the private plaine or helicopter, the sub- 
sonic jet and the supersonic transport we will derive in this chapter the 
magnitude of each error source that causes an indicated velocity error of 
1 km/h. Thereby we will first assvime that it is possible to measure the 
attitude of the gyropendulum in the navigational frame where attitude aud 
azimuth angles are zero. Attitude errors and the effect of an erroneous azi- 
muth information will be considered in a separate section. In a final section 
we will consider briefly the frequencies of the error oscillations. 
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7«1 Computational Velocity Errors of a Gyropendulum vd.th Signal Generators 
Aligned vd-th. the Navigational Freune 

The derivation of the corresponding error equations can be found in 
Appendix C, The inaccuracy terms which are the differences between the 
computed or measured terms (denoted with an asterisk) and the accurate ones, 
are indicated by a preceding (I), e.g. 


,_x-,cmd „cmd* , 

(I)Mn = Mjj + mr b^ 

,T^^u^cmd vrcmd* , 

(l)Mj. = Mj. - mr bj^ 


(7.1) 

(7.2) 


(s. eqs. (Cl) and (C2)). The initial alignment inacciuracies according to 
eqs. (C3) and (C4) are 


- ®No = 


(7.3) 


= ®Eo - ®Eo = ®e; - <=“ 


i7A) 


These values can be easily introduced into eqs. (A 38 ) and (A39). It is more 
difficult to show the effect of an erroneous angular momentum which in- 
fluences the tuning condition and the coning frequency. The computed tuning 
condition is based upon the equation (C5) 


hV(ot H*) = 


( 7 . 5 ) 


the evaluation of which may be undertaken in similar fashion to the deri- 
vation of eqs. ( 6 . 1 ) to ( 6 . 9 ) resulting in (s. eq. (C6)) 


H* = H +AH = H (1 + 
o o \ 


AS? 


1 


. 2 . + ^ 

2U)« 


r)- 


(7.6) 


This can be solved ntimerlcally when we assume that the inltiail amgular 
momentum is accurately known from laboratory tests. The eiltitude h is a 
meeisured quantity and also S2^ is based on measured Information. According 
to eq. (e?) it is 


. b » 

< ® [(u)-cosL ^ 

S cr si o 


Oi. 


.ie 




W. 


R*cos 


V • V* 

— ) AL + — tan L ] . 
h E* ^ 


( 7 . 7 ) 
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The inaccuracy within this equation which affects the computation of the 
angular momentxim is, according to eq. (C9) 

, (i)b 1 (i)v 

3(I)h V ^ 

e L rs 


w 


R 




0) 


0 _ 2 , 

R cos L 
o o 


) (I)L + 


E 


R 


tan L 




(7.8) 


The error in the computed angular momentum according to eq. (C12) is 


H 


(Dfi T.T 


(I) H = H* - H « [ 


(0 


l£ (I)h 

R 


!• 


(7.9) 


The actual coning frequency Q*^(H*) of the gyropendulum depends on the one 

hand upon the accurate values of the flight condition (h, b , L, etc.) 

and on the other upon the computed angular momentum H . Its relationship to 

c c ^ 

the theoretical frequency and the computed frequency Q is found in eq. 
(GI3) and (CI6) to be 

q'^(H*) ~ ° - (jO® (I)H/H 

o o 

a fic* s f 2(I)H (I)h 

" o I H " R 
o 0 

where can be computed similarly to eq. (6.II). 


(7.10a) 

(7.10b) 


Based upon these errors, the error of the indicated attitude of the gyro- 
pendulum could be obtained from eq. (A3^) plus (A58) and (A35) plus (A39). 

The result can be seen from eq. (C2l) plus (024) and (023) plus (025). In 
the signal flow from the indicated attitude of the spin vectors of the two 
gyropendulums to the computed ground speed some additional errors might creep 
in which should be mentioned briefly with reference to Fig. 4.3. (Remember, 
for readout in the navigational fraune: 4* = "^ = qj = 0.) 

First,, there "is the readout inaccuracy of the signal generators, which we 
denote as' (1)0- and (1)0-. The second error source is the bias error (s.eq. 
(028) and (029)) due to the difference between the computed and true values 
for the quasi static equilibrium 0^^^ , and 0^^ of the gyropendulums defined 
in eq. (4.4) and (4.5). The third and final error might occur in the compu- 
tation of the easterly ground speed from the easterly velocity with 
respect to inertial space. This error is due to erroneous altitude and lati- 
tude information (s. eq. (C4l)). 


In eq. (C32) to (C35) we have introduced some abbreviations which have the 
physical meaning of dynamic velocity error for perfect initial aQ.ignment 

(i)Av^ = 2 % * r 

o o 

(I)AV^ = 2 AVg T ^ + AM^) (7.12) 

o o 

(I)Av^g = 2 AV^g I ^ (7.15) 

o o 

where the last equation is the sum of eq. (C3^) and (C35). 

With these abbreviations we find in eq. (C36) and (G42) the following compu- 
tational velocity errors for the north and east component based on the output 
signal of either gyropendulum whose signal generators are aligned with the 
navigational frame; 

<I)V„ = (I) v;, T (I)0„ i 

- [(I)AV^ ; E^co= (I)0^*1 eosO'=(H*) t 

; [(I)AV.' iR oj® (1)0 + ] sinQ‘"(H»‘) t (7.14) 

I in o o no '' 

(l)Vg« (l)AVg + r'o)^® [cos 

o o 

* (I>®E - “o 

- [d)AvI„ i R 0)® d)0„'*’ I cos52 t 

^ X£i O O uO 

+ [d)AV^ - RghJ® ] sin^^®(S*^) * • (7.15) 

In the proposed navigational system ground speed is computed from the 
difference of the output signals of the two gyropendulums with their spin 
vectors pointing up and down. The computationeJ. velocity error has in reality 
on the right hand side of eq. (7.1^) and (7.15) » the mean error terms of both 
sensors. This dependence of the velocity error upon the mean inaccuracies has 
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especially a favorable effect upon (l)H (s. eq. (7.8) and (7.9))» But we have 
to keep in mind that the coning frequency is not exactly the same for 

both gyropendulums (s. eq. (^.3) and Table 6.l) and oscillations with oppo- 
site signs will not always cancel each other in the velocity computation. 

This might confirm us to neglect this favorable effect in the following. 

From eq. (7.11) to (7.15) it can be seen that the inaccuracies may be 
siommarized into 5 types of error sources which we will evaluate in the follo- 
wing in order to find the magnitude of each error source which yields a compu- 
tational velocity error of (l)V = 1 km/h. 

The first type of error source is the initial spin vector alignment eind 
readout error (1)0 which amounts to: 


(1)0 

TiTv 


1 


R U)® 
o o 


= 7.24 


sec 

km/h 


(7.16) 


The second type of error source is the latitude error (I)L in eq. (7.15) 
which amounts to: 


(I)L 

(iTv 


1 

R sinL 
o o 


0.0048 


deg 

km/h 


(7.17) 


at L 



This corresponds to a position error of: 


(I)P 

TiTv 


0.54 


km 

km/h 


(7.18) 


The third type of error source is the aG.titude error (l)h which at L = 45° 
is the same as the position error 


(I)h 

TiTv 




xe T 
u) cosL 


0.54 


km 

km/h 


(7.19) 


The fourth type of error source is the change of the disturbamce torques 

or the inaccuracy in the compensation of the initial disturbance torque 

^ cnid 

(I)M^ and it includes also the inaccuracy of the commemd torques (I)M 


AmVh (i)m''“‘^/H 

“(i)v - (I)v - (I)v 



0.009 


deg/h 

km/h 


( 7 . 20 ) 


0212 d 

As to the inaccuracy (l) M of the command torque, it is of interest how 
this is related to the inaccuracy in the computation of the Coriolis amd 
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centrifugal acceleration in horizontal flight, whose torque is to be 
compensated by the command torque. We find 



1_ H_ 
b, mr 




( 7 . 21 ) 


The small figures in Table 7.1a show us that the system is rather sensi- 
tive to the inaccuracy in the computation emd compensation of the Coriolis 
and centrifugal acceleration. Remember, the larger these numbers the more 
careless one might be in the con^ieasation of this effect. 


Table 7.1a Error Terms, Affecting the Operation of the Gyropendulum smd 
Causing Each a Computational Ground Speed Error of 1 km/h 



Error 

Tenn 


(i)Q 

(I)L 

(DP 

(I)h 


(i)mVh' 

<y o 

(i)m°“Vh 

(I)H/H 



Private 
Plane , 
Helicopter 


V = 250 km/h 
h = 1 km 


7.24 

0.0048 

0.54 


0.009 


Aircraft 


Subsonic 

Jet 



0.035 


7.24 

0.0048 

0.54 


0.009 


0.98 

0.024 


Supersonic 

Transport 


’ = 2,500 km/h 
. = 20 km 


7.24 

0.0048 

0.54 


0.009 


0.24 

Q.014 


The fifth type of error source is the inaccuracy of the angular* momentum 
(I)H or rotor speed control. Its effect is largest for an easterly flying 
vehicle (first term of eq. (7.12) plus second term of eq. (7.15)) 


(I)H/H 


2 (RoJ cos L + Vg) 
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which is evaluated in Table ?.1a along with the other tjrpes of errors. 

In Table 7.1b those error terms are evaluated which according to eq. (7.8) 
and (7*9) make up the computational inaccuracy of the angular momentum con- 
trol (I)H. The numerical values are derived by equating each term of eq. 
(7.8) and (7.9) with the value 2 (I)H/H^ listed in Table 7.1a. So each of 
the inaccuracies in Table 7.1b contributes a computational velocity error of 
1 km/h via the angular momentum control. 

Table 7«1b Error Terms, Affecting the Gyropendulum' s Rotor Speed Control 
and Causing Each a Computational Ground Speed Error of 1 km/h 


No. 

Error 

Term 

Dim. 

Aircraft 

Private 
Plane , 
Helicopter 

Subsonic 

Jet 

Supersonic 

Trainsport 

V = 250 km/h 
h = 1 km 

(I)H/H= 0 . 035 %''^ 

V = 900 km/h 
h = 10 km 

(I)H/H=0.024^''^ 

V = 2,500 km/h 
h = 20 km 

(I)H/H=0.014%'‘^ 

■ 

-d)h 

km 

4.460 

3.060 

1.786 

H 



86.0 

15.5 

1.87 

H 

i(I)L 

deg 

0.67 

0.26 

0.07 


-d)p 

km 

75. 

29. 

7.8 

■ 


km/h 

16.6 

11.4 

6.4 


1 ) 

s. Table 7.1a 


The values listed in Table 7.1a require indeed today's high performance 
technology; they are achievable, however. Fortunately the reqiiirements for 
the angular momentiim control are less severe. 


7.2 Attitude Errors and Req;xirements for the Azimuth Gyro 

The computational velocity errors in Section 7.1 were derived under the 
ass\imption that the signal generators of the gyropendulums are always aligned 
with the navigational frame, so errors due to inaccurate attitude and azimuth 
information were not yet included. The attitude computation is based upon the 

















output signals of the two gyropendulums and hence coupled with the velocity 
computation and its errors. From eq. (3»1l) we conclude that a velocity error 
of (l)V = 1 km/h corresponds to an attitude error of 

(I)^ _ (I) 1 _ 7 04 r7 ox'i 

TTTv “ TITv “ K ■ (7.23) 

which is not too weighty. All the velocity error contributors listed in 
Table 7.1a and b can hence be easily converted into attitude error contri- 
butors. In eqs. (^.22) and (4.23) it was shown that the vertical velocity ph 
causes an additional, attitude error which for the easterly flying SST amounts 
toA-&(ph) = 3*4 sec per ph = 1 m/sec. This is negligible for attitude control 
purposes. 

The requirements for azimuth accuracy cam be derived from a physical, in- 
sight into the system. We have seen in Chapter 3 that the gyropendulum tilts 
from the vertical in the plame normal to the motion of the vehicle. Accor- 
ding to eqs. (3*9) to (3«11) the angle is proportional to the velocity of 
the vehicle iidth respect to space, i.e., it is largest for am easterly moving 
vehicle. Since the readout taikes place in the body frame, one has to know 
true north- (the angle ) in order to evaluate this information for naviga- 
tional purposes (s. eqs. (3.20) and (3«21)). If (p is not accurately Iqiown 

(I)i|; = , (7.24) 

the system evaluates this error as a velocity in the orthogonal direction, . 
i.e., for the easterly moving vehicle as a northerly component 

(I)V^ = (I)ij; . (7.25) 

Requiring that this error shall not exceed 1 km/h, we find the admissible 
inaccuracy for (l) , which is evaluated in Table 7.1c for the three types 
of airplanes. Eqs. (5.4) and (5.5) show us that the values (l)4> are also 
applicable to the requirements for initial eJ.ignment and readout accuracy 
(1)0 of the azimuth gyro, the requirements for the drift rate and the compen- 
sation accuracy for the attitude angles (l)i&^ and (I) . 

As to the effect of gyro drift on the computational velocity error, we 
know from eq. (5.4) that this drift causes an effect about the z - (down) 
axis of = (jo'^/CjO^^ where + Vg/(RcosL) (s. eq. (5.2)) and via0jj 

causes am azimuth error of 

(1)4; = Oj^cosL = (7.26) 
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Table 7«1c Error Terms, Affecting the Operation of the Azimuth Gyro and 
Causing Each a Computational Ground Speed Error of 1 km/h 





Aircraft 

No. 

Error 

Dim. . 

Private 
Plane , 
Helicopter 

Subsonic 

Jet 

Supersonic 

Transport 


Term 


V = 250 km/h 

V = 900 km/h 

V = 2.500 km/h 




h = 1 km 

h = 10 km 

h = 20 km 

■ 

(I) ^ 
(I) 0 

min 

2.40 

1.65 

0.94 

B 

d 

0) 

deg 

h 

0.009 

0.009 

0.009 


s. eq. (A1^)). So we obtain for the azimuth gyro drift that caixses 
(I)V = 1 km/h computational velocity error ^ 


(if- _ \e (I)4> _ ± 

TTPi ~ R TiTv - R 


0.009 


km/h 


(7.27) 


This error is identical with the one found in eq. (7»20) for the gyropend- 
ulums in terms of (I)M /H. 

As to the computation accuracy requirements for the attitude angles (l)d 
and (1)^5 , we have seen in eq. (7.23) that they meet those in Table 7.1c. In 
Chapter k it had been assumed that the attitude errors caused by the vertical 
velocity have a negligible effect on the azimuth and velocity computation. 
This remains still to be proved. This attitude error is laa*gest for am east- 
erly moving vehicle (s. eqs. (4.20) and (4.21), third term on the right side) 
when AGjj(ph) according to eq. (4.22) corresponds to a pitch error A-&(ph) . The 
resulting azimuth error is (s. eq. (3.23)) 

Aip(ph) =A-9'(ph) tan L . (7.28) 

Using eqs. (7.25) j (7.28) and (4.23) we find for the resulting velocity error 


AV(ph) _ Av(ph) . A(jj(ph) 
ph ~ Aq)(phy AM ph) 


A-&( ph) 
ph 




tan L = 0.006 

m/sec 


(7.29) 


for the easterly flying SST at L = 45°. The position offset after a climb 
to a certain auLtitude h is 
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(7.30) 


(KV.„)^ tan L = 1.6?5 * 10"^ - 
i£ m 

for the same airpleine. It is in general a negligible error. If not, it can 
be easily biased out as shown in eqs. (4,2^) amd (4.25). 

7.3 Frequencies of the Error Oscillation 

Britting has shovm in [l3 ] i Chapter 2,4 that a verticaG. indicating system 
if stationary with respect to the earth, for instamce, carries out an error 
oscillation with a 24 hour period amd a Foucault modulated Schuler frequency 

(s. [13], eq. (2.31)). 

In the navigational system in question, the 24-hour error oscillation is 
generated by the azimuth gyro as shown in Chapter 5. An equation correspon- 
ding to the Foucault modulated Schuler frequency in the bracket of fl3], eq. 
(2.31) might be derived, if the difference and the sum of the gyropendulum's 
readouts are calculated in the Laplace domad.n. The common denominator of the 
combined readouts is 

(s^ +Q®^) (s^ ) = 

s + 2 s (O) +0) sin L) + (O) - u) sin L) , (7.3D 

This proves that also in the dynamics of the navigational system in question, 
no major differences are to be expected as compared with the conventional 
vertical indicating system. 
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8. INITIAL ALIGNMENT C0NSIDE3ATI0N 


In this chapter two initial alignment procedures shall be briefly outlined 
and a simplified error analysis presented. These two procedures are: 

a) The alignment based on external measurements for attitude and azimuth of 
the stationary vehicle. This allows us to draw conclusions on the actual 
stationary disturbance torques acting on the gyropendulum and to bias them 
out, 

b) The self-alignment, where based upon the initial state of the two gyro- 
pendulums and its precomputed vaLues including the disturbance torques, the 
attitude amd azimuth of the vehicle are defined. 


In these two alignment procedures, the azimuth gyro is not actively in- 
volved because it is not a north-seeking sensor. Once the direction of the 
earth axis is found, either on the basis of external measurements (item a) 
or self-aJ.ignment (item b), the azimuth gyro is slewed into this direction, 
which is then stored within the gyro. 


The basic equations for the analysis of the initial alignment are the 
equations of the accurate initiail static equilibrium about the north and east 
axis (s, eqs, (4,4) to (4,9)) 


0 . 


NoI,II 


^oI .II ~ ^0l ill 




h' lU® 

O O 


( 8 , 1 ) 


-r-r + (1) COS L_ 


EoI,II 


- ^oI.II'^'^I.II ^ 


Q 


Ol,II 


_ ^oI,II - ^ 

vr ^ _i_ _ 


xe 


cos L 


H' (0 
o o 




oI,II 


( 8 , 2 ) 


where (from eq, (6,3)) 


Q 


c 

Ol,II 



+ 



sin L 

o 


(8,3) 


They are also the equations of the actual readout angles of the x and y 
signeG. generators 


®xoI,II 


= ( 1)0 


xI,II 




- 0 


NoI,II - ®Ll,II 


(8,4) 
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(8.5) 


0 , 


yoI,II 


= ( 1)0 


yi,ii 


+ ■&, 


" 0^ T TX COS4> 

Eol ,II o 


+ 0^ sin(|J . 

NoI,II ^ o 


The asterisk denotes that these actual readout angles differ from the accu- 
rate ones in eqs. (A^5) and (A^6) by the readout inaccuracies (l)0^ 

(l)0 (which may also include errors due to vibration of the vehicle) and 
the true tilt angles 0^^ and of the spin axis with respect to the ver- 
tical, These true tilt angles again differ according to eqs. (C5) and (C4) 
from the true equilibrim in eqs. (8.1), (8. 2) by the inaccuracy terms (l)0j^Q 
and (I)0^_ . 


8.1 Initial Alignment Based on External Measurements 


If sensors for measuring the vertical (e.g, a tiltmeter for measuring 
6uid -&^) and the north direction (e.g. gyrocompass for measuring t|) ^) are avai- 
lable, one might compute from the actual readout angles in eqs. (8.4) and 

(8.5) the static equilibrium and 0^^ for both gyropendulums . The signal 
flow for this computation is shown in Fig. 8,1, Once 0^^^ and 0^,^ are known, 
they may be used directly as input signals 0^j^, and 0^^ (s, eqs, (4,4) and 

(4.5) ) in the signal flow diagram for the navigational computation (s. Fig. 
4,3) thus compensating automatically for the constant disturbeince torques 
(s. eqs. (8.1)) and (8.2)). The model which we use for the readout angles is 
similar to eqs. (a 45) and (A46) 


0 


*xol,II " ^o 
®yoI,II "" ^*o 


“ ©Lt tt cosip^^ - 0* sinijj* 


NoI,II 


Eol, II 


- 0* cos (b*" + 0^ Sind)*’ 

Eol, II ^o NoIjII ^o 


( 8 . 6 ) 

(8.7) 


where 9* » o external measurements, differing from the 

accurate ones by the sensor inaccuracies, e.g. 


(I) 9 = 9* “ 'P etc. (8.8) 

o o o 

®No ®Eo are the equilibrium angles which we want to compute; they are 
erroneous by, e.g. 


( 1)0 


No 




(8.9) 


Equating eqs. (8.6) and (8.7) with (8,4) and (8.5), we find for the error of 
this computation 
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.Bg iBltlal Allgnaent Based on External Measurements 








* l<l>®yl, II - '^^*0 I ““♦o 0 ®E0I,II (&.M 

'I>®E0l,II = '">®Ll,II - l'^>®yl,ll - 

- [d)©3^j jj ” (I)'Pq] siEUp^ . (8.11) 


In eq. (4.9) we have found that a disturbance torque expressed in 1 MERU = 
0.015 deg/h causes the gyropendulum to settle with a displacement angle of, 
12.2 s^ from the theoretical equilibrium. With respect to eqs. (8.10) and 
(8.11) we may reverse this statement. If we allow each term in these equa- 
tions only to be smaller than 12,2 s^, then each term contributes an uncer- 
tad.nty of less than 1 MERU in the final computation and compensation of the 
constant initial disturbance torque. Briefly: 


Initial misalignment 
of spin axis 


Readout 

(i)e. 

inaccuracy 

(i)©^ 

Tiltmeter 

(i)d 

inaccuracy 

(l)^^ 

With 0„ a 2.4° we 
Eo 

find from 




• ( i)m‘ 

< 12.2 sec per MERU 


( 8 . 12 ) 


ment of the gyrocompass 


(i)4j 

O 


f 

- 4,9 min per MERU 


(I)M^ 

H 


(8.13) 


Compared to the inaccTiracy for initial misalignment and readout of 7*24 sec 
contributing a velocity error of 1 km/h (s. Table 7»1a), the requirement of 
eq. (8,12) is less stringent. The achievable accuracy of sensors for measuring 
the vertical and the north direction is better than required for eq. (8,12) 
and (8.13). 
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8,2 Self Alignment 


This initial alignment procedure makes direct use of the signal flow in 
Fig, 4,5, where the input signals for a stationary vehicle are the actual 
readout angles j-j a-nd jj of the two gyropendulum's signal genera- 

tors (s, eqs. (8,4) and (8,5)) and the precomputed values for the equilibrium 
0*oi jj and jj (s, eqs. (8.I) and (8.2)). For a stationary vehicle these 

angles correspond to the angles jj and ©£<^j mentioned in Fig. 4.3. 

For perfect values of ®yoI,II ®NoI,II’ ®EoI,II Perfectly 

known azimuth 4) , the indicated initial velocity of the stationary vehicle 

must be zero, i.e., in Fig. 4.3 the differences ^©^j. - ^®yl ” ^®yll 

must vanish. This condition may be used for the self-aiignment procedure. 

Based on the inputs mentioned above, the azimuth ({j'^ entering on the left hand 

side of Fig. 4.3 first column of sin ip and cosip boxes is modified until the 

differences A0 _ - A0 and A0 ^ ~ A0 tt vanish, 
xl xll yl yll 

The angles AO^^j jj and AOy^j jj which in this chapter replace the angles 
A0^j JJ and A0yj jj are according to eqs. (4.29) and (4.3P) 

®Noi,ii + ®Eoi,n 

“ * '(’o ■ [ ^^'®NoI,II " 

■ l^^'®EoI,II ■ ^^'®EoI,Ill 
*®)oI,II “ ®)oI,II * ®EoI,II "’='*'‘0 -®Ll,II 


o ■ ®EoI,II ‘=“1', 


(8.14) 




'fo - ®EoI,II 


(8.15) 


Denoting with 6 the difference of corresponding values of the two gyropend^ 
ulums, for instaince 


6A0* = A0* - A0* 

xo xol xoll’ 


(8.16) 


we obtain for matched self ailignment condition: 
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I. 


6A0* = 0 ~ 

xo 


8 ( 1)0 - 




6(I)0j^^ ] cos (JJ 




6(1)0^^ ] sin 


+ ( I) lb 60 cos 11) 
^o Eo ^ o 


(8.17) 


and 


5 AO* = 0 
yo 


6(1)0^ - [6(I)0Eo - S(I)0t,„1 cos4) 


Eo 




.- (I)4> 80_ sinij) 
o Eo ^o 


(8.18) 


where in connection with the supposedly small error (I)tJ)^ we may approximate 


xe 


2oj cosL 

60 ~ ° 

0)^ 


(8.19) 


(s. eq. (8.2)). Multiplying eq. (8.1?) with cos(j)^ and (8,l8) with sin t|) ^ 
and subtracting both results gives for the azimuth error 


(DtJj a 


0) 


2o)^®cosL 


[6(1)0 cos4) - 6(1)0 sin4) - 6(1)0^ + 6(1)0,, ]. 

^ X o y o No ^o 


( 8 . 20 ) 


According to eq. (8.1)-, 6 (1)0^^^ is made up of the difference of the distur- 
bemce torque inaccuracies for the north axis of the two gyropendulums 


H (X) 
o o 


( 8 . 21 ) 


where again 


= <0 - < 


( 8 . 22 ) 
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is the difference between our knowledge of the disturbance torque and the 
actual disturbance torque. So we find finally 




2H' (jO^® cos L 
o o 


U) 


2af®cosL 


[6(1)0 


cos U) 

X o 


6(1)0 sintp - 6(I)0;J; 1 
y ^ o No ' 


(8.25) 


In a nonpendulous gyro, m'^/H corresponds to a gyro drift, and the first term 

on the right hand side is well known among experts. It is inherent in all 

self-aligning inertial platform systems (s. [l3]» eq. (2.45)). The terms 

within the bracket are firstly the readout errors 6(1)0 smd 6(I)Q of the 

X y 

signal generators which may also include errors due to vehicle vibration, ^d 
secondly the error 6(I)0j^^, indicating that the gyropendulums may not have 
settled to their static equilibrium about the north ^lxis (s. eq, (C3)). 

Asstiming the worst case for the alignment procedure, that all the inaccu- 
racies of the sensors are of equal magnitude but of opposite sign, we find 
for the azimuth aulignment error per 1 MEKU 


(I)4j^ 4.9 min per MEBU 


H' 

o 


(8.24) 


which corresponds to the result in eq. (8.13). As might be remembered this 
result read: the azimuth error (I)'Pq of the gyrocompass must be less than 
4.9 mj n if one wants to define the disturbance torque with sin accuracy corres- 
ponding to 1 MEBU. 

In Chapter 7, eq, (7.16) we have found that the signal generators and the 
initial, settling of the gyropendulums to their eqiiilibrium must meet the 
requirement (1)0 - 7*24 sec in order to keep the corresponding computed velo- 
city error smaller thaui. 1 km/h. Assuming the same inaccuracy during initial 
alignment, each term in eq. (8.23) nay contribute an additional azimuth ejrror 
of: 

(l)0 = 2.9 »iu per 7.24 sec readout and spin vector alignment error. (8,25) 


Once the azimuth is adjusted in the way described above, the attitude smg- 

les ,ip* and -6* are indicated in Fig. 4,3 by the sums: (A0 t + A0 tt)/ 2 and 
o o o u xoll 

(A0ypj ^^yoll^'^^* Agauji smd -&* contsd.n the inaccuracy terms, which if 
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in connection vd.th (l)4> we neglect the mean angle both sensors we 

obtain from eqs. (8,1^) and (8,15) 

“ ‘I)®x - [‘I>®No - (I)®No ! 

(I)a„ = (1)6 - - (1)6^0 ' n 




No 


(8.27) 


Again, if the readout inaccuracy and initial spin vector misalignment (l)©^^. 
and (1)0 meet the requirement of eq. ( 7 » 16 ), the corresponding attitude 

y 

error is of the same magnitude: 

(l)ip^ = (I)'9 q = 7.24 sec per 7.24 sec readout and spin vector alignment error. 

( 8 . 28 ) 


Finally an inaccuracy in our knowledge of the stationary disturbance torque 
jrrespor 
error of 


corresponding to 1 MERU causes via (^)®jjq ^ additional attitude 


(Dm" 


(Dip = (I)^ = 12.2 sec per MERU 

r\ r\ ^ • 


(8.29) 


which could have been derived directly from eq. (4,9) • 


These results show us that it is much more difficult to achieve good accu- 
racy with the self-alignment procedure compared to the alignment based on 
external measurements. 
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9 . SUI'lMARY 


An inertial navigation system is described cind analyzed based on three 
two-degree-of-freedom gyroscopic sensors mounted on a common base which is 
strapped down to the vehicle. Two of the sensors are Schuler-gyropendulums 
with their spin vectors pointing up and down. The third sensor is an azimuth 
gyro with its spin vector pointing northerly parallel to the earth axis. 

The two gyropendulums are each base motion isolated about their two input 
axes by meajis of two gimbals (one of which being the gyro case) and a gimbal 
follow-up control (s. Chapter 3, J’ig. 3.2). 

iVhen moving over the surface of the earth, the spin vector of each properly 
tuned gyropendulum will tilt with respect to the vertical in the plane normal 
to the motion with an angle proportional to their velocity with respect to 
inertial space assuming that no error oscillations are superposed on the use- 
ful signal. Looking into the direction of the motion, the arrowheads of the 
upward and downward pointing spin vectors will both tilt to the left, indi- 
cating a negative and positive angle of 7.24 sec/(]<m/h) from the vertical 
(s,. ecfs. ,(3.9) to ( 3 . 11 ) for principles and eqs. (4.35) to (4.39) for more 
details). 

Since the velocity of any point on the surface of the earth is known, the 
difference of corresponding, l*eadouts of the two gyropenduliims can be used for 
ground C.j:)«ed and position computation. The sum of corresponding readouts can 
be used' for tracking the vertical. The signal flow of this computation is 
shown in Fig. 3*5 for principles and Fig. 4.3 for more details. It is based 
on small angles for pitch and roll, but a large yaw angle which is measured 
by the azimuth gyro. 

The azimuth gyro may have a gimbal follow-up control similar to the gyro- 
pendulums (s. Fig. 3 . 3 )* It is then base-motion isolated about its two input 
axes. Unfortunately, for high latitudes the maneuverability of the vehicle is 
limited for a gimballed gyro. If latitude plus pitch angle or roll angle 
amount to 90°, gimbal-lock may occur. This is prevented by an electrostati- 
cally supported azimuth gyro where direction cosine pattern readout is used. 
Chapter 5 treats this subject. 

Error oscillations of the gyropendulums, inaccuracies of the signal genera- 
tors and the computationaul process and errors in the measurement of the azi- 
muth cause computational velocity errors. A simplified error analysis for the 
navigational system moving over a spherical earth was carried out in Chapter 7 
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v/here the accuracy requirements for the major sensitive items in the opera- 
tion of the system were defined. 

It v/as assumed that each item may contribute a computational velocity error 
of 1 Inn/h, when the system is mounted on three types of aircraft, a private 
plane or helicopter (V = 250 km/h), a subsonic jet (V = 900 km/h) and a 
supersonic transport (SST, V = 2,50Q-km/h). The major sensitive items in the 
operation of the system are: 

a) The gyropendulum must be properly tuned in order to prevent horizontal 
and vertical accelerations from starting oscillations with approximately the 
Schuler period of 84 minutes (s. Chapter 4). The tuning condition (s. eqs. 
(4.10) and (4.3)) is a function of easterly velocity, latitude and vertical 
Coriolis and centrifugal accelerations and must be matched by a proper con- 
trol of the angular momentiim (s. Fig. 6.1). The change in angular momenttim 

is approximately 6 % for an easterly flying SST at 40° to 50° latitude (s. eq. 
(6.9) and Table 6.1). An inaccuracy of 0.0l4 % in this angular momentum con- 
trol causes a velocity error of 1 lon/h on the SST (s. Table 7*1a). 

b) The Coriolis and centrifugal accelerations due to horizontal velocity 
must be properly compensated by means of command torques, applied to the 
gyropendulum (s. eqs. (4.12), (4. Ip) and Fig. 6.1). Expressing the command 
torque as M° /H, an inaccuracy of 0.C09 deg/h causes a computational velo- 
city error of 1 lim/h (s. Table 7»1a). This means that the Coriolis and centri- 
fugal acceleration must be compensated to an accuracy of 0.24 ;o on the SST. 

c) The gyropendulums must have settled initially to their static equili- 
brium, v;hich in the absence of disturbance torques is 0° from the vertical 
about the north axis and approximately + 2.4° (+, because spin vectors are 
pointing up and down) about the east axis. The latter value is due to the 
equilibrium of the gyroscopic torque caused by the horizontal component of 
earth rate and the pendulous torque (s. eqs. (4.4) to (4.9)). An initial off- 
set angle of 7.24 sec. between the spin vectors of the gyropendulums and 
their equilibrium causes a computational velocity error of 1 km/h (s. eq. 
(7.16) and Table 7.1a). Constant disturbance torques, causing constant gyro 
drift in a non-pendulous gyro, cause a constant bias from this equilibrium 

in the gyropendulum. This may be compensated if knovm beforehand or measured 
during initial alignment (s. Chapter 8). 

The compensation error of the disturbance torque M°’/H must be smaller than 
0.009 deg/h in order to keep the resulting computational error below 1 lun/h 
(s. eq. (7.20) and Table 7.1a). The same figure applies to the change of the 
dusturbsmce torques during flight, be it deterministic or random in nature. 
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d) The velocity computation is based upon the output signals of the four 
signal generators in the two TDF gyropendulums . A readout inaccuracy of 
7.24 sec corresponds to a velocity inaccuracy of 1 km/h (s. eq. (7.16) and 
Table 7.1a) 

e) For the computational process to obtain ground speed from the output sig- 
nals of the gyropendulums, the azimuth with respect to north must be known. 

On an easterly flying vehicle, for instance, any azimuth error will be inter- 
preted by the system as a northerly component. The azimuth accuracy is effec- 
ted by initial alignment accuracy and gyro drift, both causing the gyro on a 
stationary vehicle to oscillate with a 24 h mode with respect to the earth. 

It is effected furthermore by the readout accuracy and compensation for the 
vehicle's attitude. In order to keep the computational velocity error smaller 
than 1 Ina/h, the accuracy requirement for the azimuth (alignment, readout, 
attitude compensation) must be knov.n within 0.94 min on an SST (s. eq. (7.26) 
and Table 7.1c). The corresponding requirement for gyro drift is 0.009 deg/h. 

f) Since, in the computational process for ground speed, the earth's velo- 
city has to be subtracted, latitude and altitude of the vehicle must be icnown. 
The former is an output signal of the system itself (s. Fig. 4.3), the latter 
must be measured externally. For a computational velocity error of 1 km/h 
latitude and altitude, i.e. position and altitude must be known within 0.54 Ion 
(s. eqs. (7 .i 3), (7.19) and Table 7.1a). 

This study includes some design considerations in Chapter 6 and initial 
alignment considerations in Chapter 8. 


Dr.-Ing. Bernhard STIELEB, Institut fiir Flugfiihrung 

der DFVLR, 

33 Braunschweig-Flughaf en 
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Appendix A 


Derivation of the Performance Equations for a Gyropendultun with Vertical Spin 
Axis in a Dynamic Environment 

A1 Readout Angles in the Navigational Frame 

The model of a pendulous two-degree-of-freedom gyro with vertical spin axis, 
used by Schuler in [l] is 'shown in Fig. A1 . We will assume that the gimbal 



02 

A1 Modal of the Pendulous Two»Dsgraa»of»F>>aedeD Oyz^soope with 
Vertical Spin Axis 

A spherioal gimbal bearing 

rotor bearings 

S center of gravity 

r lever arm of pendolosity 

bearing A in this model is nearly frictionless, so that the rotational motion 
of the base on which this bearing is mounted, exerts only a small disturbance 
torque on the gyro. The attitude of the gyro's spin vector with respect to 
the vertical is affected mainly by the pendulous torques due to linear accel- 
erations. Let us assume for the present that we may read this attitude in 
the navigational fraune (subscript "n") about the north and east axis as shown 
in Fig. A2. 
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In the following we will derive the performeince equation of the gyropendulum 
in terms of these readout angles 0jg and 0g. 


Applying Newton's Law 


M = ^ 
dt 


dE 

dt 


+ U)^^ * 


(A1) 


S 


to the gyropendulum, we assume that in the total angular momentum H , the 
spin angular momentiun H is the dominaint part and constant 


H a H = const. 


(A2) 


So we may rewrite eq. (Al) in the gyro frame (subscript "g") as 


X H , 

S g 


(A3) 


where according to Fig. A2 

.0 


H = H 
g 


(?)• 


-♦iff 

The angular velocity uj “ may be split up into two components 


(A4) 


-*ig -^n -*ng 
U) ^ = (Ji + Ui , 


(A5) 


where the latter component is the angular velocity vector of the gyro with 
respect to the navigational frame (s. Fig. A2, 0„ and0' are small angles) 

N X< 


CO 


Jig - 



(A6) 


and the former one is the angular velocity vector of the navigational frame 
with respect to inertial space. In the navigational frame, this vector is 



By means of the transformation matrix 



(A7) 


(A8) 
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n 


may be transformed into the gyro freime 



C 



in 

n 


(A9) 


Let us now introduce eqs. (A4) to (A9) into eq. (A3) and assume that we have 
two gyropendulums, the spin vector of the first one (upper sign) pointing up 
and that of the second one (lower sign) pointing down. The basic performance 
equation for these sensors is 



(A10) 


In eq. (A7) let us write the rate of change of latitude and longitude in . 
terms of the northerly aind easterly ground speed on an earth assumed to be 
spherical 

pL = V^/R (A11) 

pi =. Vj^/(R cos L), (A12) 

where 

R = R + h (A15) 

o 

is the magnitude of the radius vector from the center of the eairth to the 
gyropendulum. The easterly velocity with respect to inertial space is: 

V = Ru)^^ cos L + V . (A14) 


The rate of change of celestial longitude is defined as 
pX = U)^® + pi = cos L). 

We introduce this and eqs. (A1l) and (A14) into eq. (A?) and obtain 
^p X cos L \ 


(A15) 


0 ) 


in 




.-p X sin 



-pL = -Vjj/R 


-(V.^R) tan L . 

Xi!i - 


(A16) 


78 


Let us now draw our attention to the torque M acting oh the gyropendulum. We 

S 

may write 

iT = mr* X r* + TT + TT (A17) 
g g S g S 

where is the disturbance torque which generates the gyro drift in a non- 
pendulous gyro, is the command torque generated by special tbrquers 

(s. Fig. 5»2), m is the mass of the pendulosity, 



(Al8) 


is the radius vector from the gyropendulum' s center of support to the center 
of gravity and f is the nongravitational specific force which the gyropend- 
ulum exerts on the gimbal bearing A in Fig. A1 . This specific force vector is 

f** = "g" - a^^® (A19a) 

where G is the gravitational field intensity vector (mass attraction) and”a^® 
is the acceleration vector of the gyropendulum with respect to inertial space. 

Since it is more convenient to calculate in terms of the gravity vector g, 
defined as 

^*=■^-0)^® X (ui^® x'r), 

we obtain for eq, (A19a) after a few additional calculations 


f“ = T- V -~T 


(A19b) 


where a is the acceleration vector with respect to an earth-fixed point and V 
is the groiind speed vector. In eq» 12.5 we find an expression for the 

specific force f^ in the navigational fraune. According to [l3], eq. (2.1), 
this is the negative of the force defined in eq. (A19)» Assuming a spherical 
earth (in [l3ls © = 0, r^ = r^^ = R) and zero deflection of the vertical 
(in [13] J ^ = 1) = 0) [l3]» eq* (2.3) reads with R = ph 

y^-R p^L - ^ r[ (pX)^ _ U)^® ^]sin2L - 2ph; pL 
-R cosL p X + 2R sinL pL pX - 2ph pX cosL 

\g + p^h - R [(pX)^ - cos^ L - R(pL)^/ (A20a) 


f = 
n 
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Collecting terms, we may derive the following equation 
/- P\ - ph Vj^K- [ V^g/R - R',(a)^®cosL)^l tan L 


f = 
n 


- tanL)V^g/R 

g + p\ - V^/R - V^|/R + R (0)^®cos3l)^ 


(A20b) 


which may be transformed into the gyro frame as was eq. (A9). 

The Coriolis and centrifugal acceleration shall comprise two different 
terms b and c, where b refers to a horizontally flying vehicle 


^ - P ((jj^®cosL)^ltan L, 


b = vJ/R + V.J/R - R (u^®cosL)^ 
V N' lE' 


(A21a) 

(A22a) 

(A23a) 


and c is the additional Coriolis acceleration caused by the vertical velocity 
ph 


®N ~ P^ 


=E = W® 


0 = 0 . 

V 


(A21b) 

(A22b) 

(A23b) 


So we can rewrite eq. (A20b) in the following form 


/ 


N 




E 

g + ^ S + 


pV^E " “ °E 


p^h _ b. 


(A24) 


We transform this into the gyro frame similarly to eq. (A9) and put the re- 
sult into eq. (A17) in order to find the torque acting on the gyropendulum 


P^iE + ^E ®E - + P H - b^) 


„cmd -,d 

+ % 


”g ' ” \- - "n - 

0 0 
Inserting this expression auid eq. (A16) into (AlO),. we obtain 


(A25) 
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L 



'mr(g + - b^) i H V 


iE^ 


i Hp 


tan L + Hp 


mr(g + p^h - b^) - H L /\0 



■mr (PV.J, ^ bj. ^ Cj.) ; HV^/E + ^ < 

-mr (pV^ + bp + Cp) * HV^p/E + 


(A26) 


Before taking the Laplace transform of eq. (A26), we have to specify the 
coefficients which are time- dependent in general , 6specially. those of the 
matrix on the left-hand side. The frequency ,found by the evaluation of its 
determinant and describing the undamped coning motion of the gyropendulum 


2^ ( -| + 2^ - —) 1 cos L + tan L (A2?) 

H g g R 

2 

depends on the vertical acceleration p h which has zero mean indeed for a 
horizontally flying vehicle, but may have short time peaks of up to 1 g or 
more. In order to estimate its effect upon the motion of the gyropendxilvun we 
assume that this vertical acceleration is an impulse function . which occurs 
at the time T 

p^h = ph 6(t - T ) (A28) 

whose product with the current angle (t) and 0g (t) becomes 

mr p^h 0JJ g(t) = mr ph 6 ( t - T ) 0^^ * (A29) 

Since 0 (t) is a constant angle, this enters the right hand side of eq. (A26), 
leaving the slowly varying frequency on the left hemd side 

fiC = (>| 1) 1 (o)^® cos L. + tan L . (A30) 

n g « 


For the singles 0^j and 0^. we will introduce a quasi static eqtiilibrium 0^^^ , 
0„^, and a time dependent dynamic part A0 and A0„ : 

rfl r« £1 


N N1 N 


Qe = ®E1 - ^®E 


(A31) 

(A32) 


For all the other variables on the right hand side of eq, (A26), we introduce 
a steady state term (Subscript "Q.") referring to the initial value smd a 
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displacement term like, for instamce: 


V = V + AV . 
o 


(A55) 


So we obtain for the quasi static equilibrium of eq. (A26) for zero command 
torque 


© 


M? 

1^0 


Mgp/H + 0)"'^ cos 

■ 


0. 


El 


(A34) 


(A35) 


leaving as equation for the displacement angle from this quasi static equi- 
librium 


Q 


; p\ /A© 


■- p fi'"/ V A© 


-^ph6(t-x) 




©E 


mr 

H 


‘P ♦ '’e °E> + ^ \ 

- |£ (p AVj, * * c„) J AV.J./R t Oq"^ * Am^)/h/. 


(A36) 


Let us now specify the coefficients and independent variables. We assume that 
R, COB L and sin L are constant. The time dependence of AV^^ and is speci 

fied only insofar as we assume that for t = o the velocity does not change 
discontinuously such small terms as the Coriolis amd 

centrifugal acceleration b- in horizontal flight, the additional disturbance 
torque Am^ and the command torque, we simply assume step functions. With the 
Laplace transform of the additional Coriolis acceleration c„ and c„ due to 

iM - . ill 

the vertical velocity (s. eqs. (A21b) and (A22b)), one has to keep in mind 
that with the assumption in eq. (A28), the Laplape transform of the vertical 
velocity isX|phJ = phe”®^/s. Since in eqs. (A21b) and (A22b) V^^/R emd 

will be assumed to be constant at the time t - T , the Lanlace trajisform of c„ 

f 1 St > N 

say isX{cj^j = c^e"® /s. 

Denoting with A the determinant of the left-hand matrix of eq. (A36) 

A= + . (A37) 

we obtain the following equations for the displacement angles from the quasi 
static equilibriiim in eqs. (A3^) and (A35) 
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A0 


A0 


A2 


1 r - mr / 2 nc H y ,mr lx 

N “ aI"^ H ■*■ mr ^ " B^ ® 

n° r /, ^iE , -ETx «Tcmd A^dl 

1 I [mr (b^ + / Ph e“®^) - - AM^ ] 

^ ^'"[^®Eo - ^ P^] 

•" ®[^®No ^ ] 1 

1 r - mr / 2 ^c H ^ /mr ^c 1^ _ 

E = 4 I* r “ 3 Tb> '"'iE - <H- ” - R> ^ ^’'h 

- S [- '^N - r P*- * < ] 

* 5 ["“■ '‘e ♦ ^ J’*' ♦ “n”"^ + '^4 ] 

- ”°[^®No rP“] 

*s[4e^„l9,(x)e-^|Eph]]. 

Beadout Angles in the Body Frame 


(A38) 


(A59) 


The eqs. (A 38 ) and (A39) transformed into the time domain and added to eqs. 
(A3^) and (A35) describe the attitude 0^^ and 0^ of the gyropendulum's spin 
vector with respect to the vertical about the north and east axis of the navi 
gational frame. If, as proposed in Fig. 3*2, the actual gyro readout taikes 
place in the body frame, we have to derive the corresponding equations for 0^ 
and© , the readout angles along the vehicle's longitudinal axis and pitch 

y 

axis. They may be obtained in the following way: 

The unit-spin vector l“ = g, given in the navigational frame, may be trans- 
formed into the body frame via the rotation (p about the vertical or yaw axis 
( (|> = 0 for northerly flying vehicle), the rotation about the pitch axis 
( •& > 0 for tail down), and the rotation <p about the longitudinal axis (<p > 0 
for right wing down) . So we obtain 





C 

ng 



(A40) 


where 
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'cos (p sin 4> 

cos 

;dcos 4> + <4> sin 4> -M>cos4> +-&sin(jj 


^b» ~ 


(A^l) 


(the attitude angles •& and <i) are assumed to be small); C can be taken from 

ng 

eq. (A8), and 



(A42) 


So we find 


1® ~ 


(-•& + 0J, cos4> - 0^ sin4>\ 
ip - 0j^ cos 4> - 0£ sin (jji 


(A^3) 


The unit-spin vector in the body frame may also be described in terms of the 

readout angles 0 and 0 , i.e, the angles between the gyro frame and the body 
X y 

frame (positive when the body moves' positively around the x and y axes). 

Since the attitude angles ^ and ip of the vehicle with respect to the navi- 
gational frame are small, as well as the attitude of the spin vector with 
respect to the navigational frame, also 0^ and 0^ are small, which are com- 
posed of both kinds of angles. From comparison with Fig. A2 and eq. (A8) 

(replace by © , ©„ by 0 , the n-frame by the g-frame and the g-frame by 
IM X ii# y 

the b-frame), we obtain 


« I 0. 



- ©„ 


0 


X 


(A44) 


Eqs . ( A43) "and (A44) describe the saime spin vector, so we may deduce: 


Og siniiJ 
©jj sintp 


(A45) 

(A^6) 


A3 Application of the Command Torques in the Free Azimuth Fraime 

Command torques which have to be applied to the gyropendulum were denoted 
in the previous sections as sind as torques about the north and 
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east axis of the navigational frame. The magnitude of these torques is defined 
in Chapter 4, eq. (^.12) and (4,13). 

In the proposed design for the gyropendulum , the gyro case including the 
gyro torquer is base-motion isolated about the input axes but not about the 
spin axis of the gyro. So the reference frame in which the torque is applied 
to the gyropendulum coincides approximately with the free azimuth frame which 
is rotated relative to the navigational fraune through the azimuth angle (j) . In 
reality there is a small difference caused by the angles and 0^. which we 
will neglect, since for the small angles 0^^ or 0^. (< 5 ) horizontal command 
torques are only affected by cos 0^^ or cos 0^, which can be assumed to be 
unity, especially since the command torques for compensating Coriolis and 
centrifugal accelerations are small. So we write: 



C 

an 


M 


cmd 

n 


where 

( cos4j sin4i 0 \ ' . 

-sin4> cos 4) 0 I 

0 0 ' 1 j 

(s. eq. (A4l) for ■& = ip = O) . - 

For the components of the gyro torquer we obted.n 


„cmd , „cmd ■ . , 

~ ^ costp + Mg sin(|) 

,.cmd «jfCffld ill rnrCmd « i 

My = Mg cos4) - Mj^ sintp . 


(A47) 


(A48) 


(A49) 


(A50) 
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jendix B 


Derivation of the Performaince Equations of a Two-Detrree-of-Freedom Azimuth 


Gyro vd.th its Spin Axis Pointing Northerly. Parallel to the Esu?th Axis 


B1 Readout Angles in the Equatorial Tangent Frame 

V/e define the equatorial tangent coordinate frame (subscript "t”) as a 
f ramie with the x or N'-axis pointing northerly pairallel to the earth axis, 
the y or E-axis pointing horizontal to the east and the z or D-axis pointing 
downward tilted from the vertical by an angle equail to the latitude angle 
(s. Fig. Bl). 



The basic equation for the performance of a TDF azimuth gyro with its spin 
axis pointing northerly pairallel to the earth axis may be obtained from eq. 
(AIO), the corresponding equations for a TDF gyro with vertical spin aixis, 
by a cyclic permutation of the subscripts as follows:: 

N-*E 

E-*-D. 

Keeping in mind that the azimuth gyro is non-pendulous (r = O), we obtaiin 
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it 

where the only input angular rate U) is 
0)^^ = U)^® + Vg/(R cos L). 


(B1) 


(B2) 


Assuming first that the only torque acting on the gyro is the constamt di- 
sturbance torque generating the gyro drift 



(B3) 


second that the gyro has an initial misalignment about the east and down axes 
and 0^^ and finally that the easterly velocity V^, is constant, we obtain 
in the Laplace domain: 


0., = 


1 /to 

A ^ s 


d 


it 


®Do * “e ♦ = ®Eo> 


it 

r\ 1 / to ,.d it n+ d r»+ N 

®D = S — “e - “ ®Eo * = ®Do' 


where 

A= s^+ (U)^^ 


(B^) 

(B5) 

(B6) 


The signal flow diagram for this gyro is very similar to the one for the gyro- 
penduliim in Fig. 4,1 except the gyro is nonpendulous (r = O) and the sub- 
scripts are different (as mentioned above). 


B2 Readout Angles in the Body Frame for a Gimballed Gyro 


As mentioned in Chapter 3» the azimuth gyro is base-motion isolated about 

its input axis by means of one gimbal between the gyro case and the vehicle 

(s. Fig. 3«3)» the gyro case serving as second gimbeil. In order to derive the 

readout equations of the two signal generators for 0 and 0 in the. body frame 

y ^ 

we proceed an in Appendix A, Section A2, i.e., we first describe the unit spin 
vector 1 in the body frame in terms of the angles 0^, 0^^, L (latitude) and(|) , 




•d, tp 



Cv C . . C.. 
bn nt tg 



(B7) 
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where C. is the trams formation matrix between the gyro frame and the equa- 
torial tangent frame (rotation through 3^, 0^^) 


^tg “ ^ ®D 


■ ®D ®E 
1 

0 


(B8) 


is the transformation matrix between the equatorieG. temgent frame amd the 
navigational frame (rotation through L about east axis) 



/ cosL 

0 

sinL 

c - ^ 

= [ 0 

1 

0 

nt 

I 




\-sinL 

0 

cosL 

and 

C, from 
bn 

eq. 

(A41). 


) 


(39) 


So we obtain 

( cos 4) cos L + ( •& - 0g cos 4 j ) sin L + 0^^ sin 

-sin cos L - ( - 0g sincp ) sin 1+0^^ cos ll> 

-sin L + (^ cos <|) + 9 sin (J) - 0^) cos L ji (BIO) 

We now derive the same unit spin vector in terms of the readout emgles 0 

y 

and 0 which are positive when the gimbal rotates with respect to the gyro 
z 

about the y-axis and the vehicle rotates with respect to the gimbal about the 
z-axis. This results in 


''cos 0 cos 0 sin 0 

y z 2 

= {-cos 0 sin 0 cos 0 

og \ y z 2 

sin 0 0 


-sin © cos 0 
y z 

sin 0 sin 0 
y z 

cos 0 


(B11) 


amd 


'cos © cos 0 
y z 

ir = I -cos 0 sin 0 
b 1 y z 

sin 0 


(B12) 


By equating eqs. (BIO) and (B12) we fiid 


cos 0y sin 0^ = sintj> cos L + ip sin L - 0^ sin(j) sin L - ©^ cos (p 


(B13) 
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and 


sin 0 

y 


sin L + ( -^cos (J) + ip sin ip “©g) cos L. 


Using the formula 


(Blif) 


arc sin (u^ 1 - i v - u^) = arc sin u - arc sin v 
and putting 


(B15) 


u = - sin L 

V = <& cos (p + vp sin Jj) - , 

we may approximate eq. (B14) with v considered small as 

0y ~ - L + -& cos 4> + sin - 0g« (Bl8) 

This angle might very well become 90° at high latitudes for which eq. (B13) 
deteriorates, indicating that ^'gimbal lock" occurs. For situations when this 
is not the case, v;e may put eq. (Bl8) into (B15) s-nd solve for sin 4) . In 
connection with small quantities, we approximate sin © » sin 4> and find in 

a way similar to eqs. (B15) to (BI8) 


(BI6) 

(BI7) 


0 _ 


cp - ( -& sin 4> “ cos 4* ) tan L + 


Qp 

cosL 


(B19) 


The next section deals with the readout equations of a gyroscope, where 
"gimbal lock" does not occur. 


B3 Readout Angles in the . Body Frame for, an Electrostatic Gyro (ESG) with 
Direction Cosine Pattern Readout 

No gimbals are required for the operation of an electrostatic gyro with 
direction cosine pattern readout (as described in [12]) which removes the 
danger of gimbal lock. 

The performeince equations in section B1 might be very well applied to the 
ESG, also eq. (BlO) for describing the attitude of the unit vector of the 
angular momentum with respect to the body freune in terms of the displacement 
amgles 0^ and 0^^, the latitude L, the attitude -9 and jp ajad the azimuth 4^ . 
Only the readout amgles which we denote with 0^ and 0^ are different; they 
are shown in Pig. B2. Instead of eq. (B12) we may write now 
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which, may be equated with eq. (B10). Similar to eq. (Bl8) we find from the 
z-component 

0y ~ 90° + L - ■& cos <i> “ *P sin 4> + 0^ ; (B2 

from the y-component we obtain 

0^ = arc cos [- sin (}) cos L - ( q>- 0^ sin (|) ) sin L + 0^^ cos (!>], (B2, 

which shows that at the north pole only ij) camnot be measured l^y 0^ . 



Fig, B2 Readout Angles of an Electrostatic Gyro 



Appendix C 

Derivation of the Error Equation of the Gyropendulvun Within the Navigational . 


jerame 


In the following we will denote all computed values by an asterisk. They 
differ from the accurate ones by the inaccuracy, e.g., (s. eqs. (^. 12 ) and 
(^.13)) 


cmd 

„cmd* , 

N 

= + mr b. 

cmd 

cmd’^ 

E 

= Mp, - mr b. 


(Cl) 

CC2) 


The initial alignment inaccuracy of the spin vector from the equilibriuin in 
eqs. (4.4) and (4.5) is 


^^^®No ^ ®No " ®No ■ ®No " 

^^^®Eo = - ®Eo = ^o>/«o • 


(C5) 

(C4) 


These inaccuracy terms may be put directly into eqs. (A 38 ) and (A59)» It is 
a little more difficult to derive the equation for the inaccuracy of the 
tuning condition. The computed timing condition is derived from the equations 


H' 


mr R* 




(C5) 


which, when evaluated similarly to eqs. ( 6 .I) to ( 6 . 9 ), gives for the computed 

.c* 


angular momentum H 


H* = H + AH* = H (1 + 
o o 


Po h* , 


where we assume the zero- values to be accurately known. The termAfi 
known quantities; An^*fg is computed from eq. (6.7) 

• n 


h* are 


H 


(C6) 


and 


1ft * b* T* V* 

a - i — [(o,^®cosL + — ) L* + -^ tan L ] . (C?) 

% Q * ft • 


0) 


g 


(b 


s 


H cos 




This computed frequency change for constemt angular momentvun, which in our 
aneilysis is the baisis for the computation of the atngular momentum H* , differs 
from the exact value by the inaccuracy 





(c8) 
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where in similar fashion to eq. (C 7 ), 


d)Q" 


(I)b. 


ID 


Ho^ -2(l)h _ V + 1_ f(^,i 

— ~ -R g 


xe 


CD' 


cosL + 
o 


2 

Rcos L 


-) (I)L + 


(I)V (C 9 ) 

— - ■■■ - tan L J. 

R O'* 


V/ith the known computed angular momentum we may tackle the computation 
of the tuning condition in the first parenthetical expression on the right 
hand side of eqs. (A38) and (A 39 ) 


mrR 


C/TT*' 




mr R^^d - (I)h/R*) 


(1 + (I)h/R ). 


(CIO) 


C ' 

The relationship between the computed coning frequency fl and the actual 
c Jic 

frequency (H ) which is based on the erroneous angular momentum is still 

C 3^^ C ^ A 

unknown. In order to compute (H*) in terms of we write for H* 


H = H + (I)H, 


(C 11 ) 


where similarly to eq. (6.9), 
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(C 12 ) 


So we obtain, with eq. (A30), the actual frequency 
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The accurate frequency, written in terms of the computed one and the inaccu- 

1 

racy term is 


- (I)fi°, 


(Clif) 


where (Dfi'" is similar to eq. (6.II) 


(l)Q‘ 




(I)H (I)h 

H R 

o o 


)• 


(CI5) 


So we find finally 

n=* = .n' (Dh 
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(CI6) 


and using eq. (CIO) for the tuning term in the first parenthesis of eqs. 
(A38) and (A39), this becomes 
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b 




(C17) 


Similarly, the second parenthesis on the right hand side of eqs. (A38) and 
(A39) becomes 


fin'CH*) - 1 = - 3 ^!j“(h‘) . 


R 




(CI8) 


Before we solve eqs. (A38) and (A39) in the time domain, including the error 
equation from eqs. (Cl) to (CA) and (C17) and (C18), we have to specify the 
transient of AV„ and AV.„. For the derivation in Chapter 4 these transients 
did not have to be specified, except for continuity (AV^^^ = AV^^^ = 0) 
because for accurate tuning, the first parenthesis, of the right hand side of 
eqs. (A38) and (A39) cancels with the determinant A and the second parenthe- 
sis becomes zero. This meant that the accurately tuned gyropendulum follows 
any transient of AV without time lag or oscillation. That does not hold any 
longer for inaccurate tuning and the transients have to be specified for 
evaluating these equations in the time domain. VJe will assume that AV . and 
Av^g are step functions but with no discontinuity at t = 0. Mathematically 
this is verified by writing for instance 


AV, 


N 


" s(l + Ts) 


(C19) 


which in the time domain is 

AVjj (t) = Vj^^ (1 - e"''^/^). (C20) 

-t/T 

V.^e assume that T is so small that for the time of observation e is negli- 

gible. This matches the physical world much better than the accurate step 
functions. 

So we find from eqs. (A38) and (A39) for the indicated displacement angles 
about the north and east axis with respect to the equilibrium in eqs. (A3^) 
and (A35) 
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(C22) 


As to the last bracket on the right hand side of eqs. (C2l) and (C22), it can 
be shovm to be negligible even if 0„(x) and 0„(x) are not exactly proper- 

rJ £j 

tional to + Vj^^/(RQ^(H*) ) and - V^^/CRQ^CH*) respectively, as had been found 
in eqs. (^.l8) and (^.19) for an accurately working gyropendulum. Assuming 
that from either amplitude within the bracket, an inaccuracy angle (l)0 
remained, the vertical acceleration would cause the gyropendulum to oscillate 
with an amplitude of p ph (l)0. In terms of indicated velocity error this 
would be 

(l)V(ph) = ph * (I) 0 . (C23) 


9 ^ 


Since ph is in the order of meters per second and (1)0 in the order of arc 
seconds, this can really be neglected. 

In Chapter 4 we have also found that the third term on the right hand side 
of eqs. (C2l) and (C22), i.e.», the bias due to vertical velocity, can be 
left out. It drops out in any case in the computation of the difference of 
corresponding readouts of the two gyropendulums, i.e., it does not directly 
affect the velocity computation. Its effect on the attitude computation, 
which is based on the siim of corresponding readouts, is negligible (s. eq. 

(4.23)). 

The total indicated angle of the gyropendulum about the north and east axes 
(s. eqs. (A 31 ) and (A32)) comprises firstly the quasi static equilibrium 
(s. eqs. (A34) and (A35)) where the denominator is now (H ) 



(C24) 


+ C0^®cosL^ 
fi°(H*) 


H CO® 
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, 'i ig 

Ho - oj cosLq 
CO® ^ Q‘"(H*) * 


It comprises secondly the dynamic angles A0j^ and A0g from eqs. (C2l) ajid 
(C22) and thirdly the readout inaccuracies (l)0j^ ^nd (1)0^. So we find for 
the total indicated angles about the north and east axes 

0* = 0 + ( I )0 = 0 + A0 + ( I )0 ( C 


E ~ E 


E “ El 


In Fig. 4,3 it is shown that the biases of the readout angles, which accor- 
ding to eqs. (4.4) and (4.5) are the angles due to initial horizontal earth 
rate and disturbance torques, are subtracted from the readout angles. This 
may cause an additional error for slightly erroneous biases which, we denote 


( 1 ) 0 . 




H*S^°(H*) H CO® 
o o 
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(C29) 
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Multiplying the readout angles of one gyropendulum by ~ gives the computed 
velocity. In connection ’.vith small quantities, we may approximate in eqs. 
(C21) and (G22) 

TT* • ■ 

» R (C30) 

mr Q ^ ( H * ) ° 


(s. eq. (Cl8)) and 

(C31a) 

o 

H*a H . (C31b) 

o 


Before we write down the computational velocity errors, let us introduce 

_ jj 

some abbreviations into eqs. (C2l) and (C22) multiplied with + — . First we 

^ . mr 

split up into its two components + AV l^dth 


AV^® = RCJ-® (cosL - cosL ) (C32) 

o 

(s. eq. (A14)). Then we define as "dynamic velocity errors for perfect initial 
alignment" 
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(C53) 


(C34) 


(C35) 


So we finally find an expression for the computational error of the velocity 
with respect to inertial space. (The formal steps for the north component . 

say are as follows: take eq. (C22), add the readout inaccuracy (l)0,, , sub- 

- H* ^ 

tract eq. (C28),- multiply by + — and take into account eq, (C30) to (C35)) 
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In this latter equation, another, error source has to be included, the error 

occuring in the computation of the easterly ground speed from the easterly 

velocity Avf with resnect to inertial space. Fig. 4.3 shows us hov/ this is 

accomplished. First the total easterly velocity V. with respect to inertial 

is* i 6 

space is computed by adding to the computed velocity = R W cosL^ 

of a. point above the starting point with the seune radius as the vehicle: 


tA©* Atr* le T atA 

V = V + AV.„ = R 0 ) cosL + AV.„. 
xE o xE o xE 


From this equation we may derive the error 

(I)V.„ = (l)h(/®cosL. + (I)AV.„. 
xii o xE 


(C38) 


(C39) 


Then in Fig. 4.3, the computed velocity R u) cosL is subtracted in order to 


obtain the easterly ground speed V^. Since 
L*" = L + (I)L, 


(C40) 


we obtain 
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sin L*’(I)L. 


(C41) 


= AV?-, + S* CO^® (cosL - cosL*) 
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oo we find the computational error of the easterly ground speed to be 
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